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LIMIT THEOREMS FOR POINT PROCESSES UNDER GEOMETRIC 
CONSTRAINTS (AND TOPOLOGICAL CRACKLE) 

By Takashi Owada t and Robert J. Adler*^ 

Technion - Israel Institute of Technology 

We study the asymptotic nature of geometric structures formed from a 
point cloud of observations of (generally heavy tailed) distributions in a Eu¬ 
clidean space of dimension greater than one. A typical example is given by 
the Betti numbers of Cech complexes built over the cloud. The structure of 
dependence and sparcity (away from the origin) generated by these distribu¬ 
tions leads to limit laws expressible via non-homogeneous, random. Poisson 
measures. The parametrisation of the limits depends on both the tail decay 
rate of the observations and the particular geometric constraint being consid¬ 
ered. 

The main theorems of the paper generate a new class of results in the well 
established theory of extreme values, while their applications are of signifi¬ 
cance for the fledgling area of rigorous results in topological data analysis. In 
particular, they provide a broad theory for the empirically well-known phe¬ 
nomenon of homological ‘crackle’; the continued presence of spurious ho¬ 
mology in samples of topological structures, despite increased sample size. 


1. Introduction. The main results in this paper lie in two seemingly unrelated 
areas, those of classical Extreme Value Theory (EVT), and the fledgling area of 
rigorous results in Topological Data Analysis (TDA). 

The hulk of the paper, including its main theorems, are in the domain of EVT, 
with many of the proofs coming from Geometric Graph Theory. However, the con¬ 
sequences for TDA, which will not appear in detail until the penultimate section 
of the paper, actually provided our initial motivation, and so we shall start with a 
brief description of this aspect of the paper. 

Many problems in TDA start with a ‘point cloud’, a collection X = {x \,..., x n } 
of points in M rf , from which more complex sets are constructed. Two simple exam¬ 
ples are the simple union of balls 

n 

U(X,r) = |J B(x k -r), 

k=i 
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where B(x: r) is a closed ball of radius r about the point x, and the Cech complex, 

C(X,r). 

DEFINITION 1.1. Let X be a collection of points in M. d and r be a positive 
number. Then the Cech complex C(X, r ) is defined as follows. 

1. The 0-simplices are the points in X. 

2. A p-simplex o = [xi 0 ,... ,Xi ] belongs to C(X,r ) whenever a family of 
closed balls { B(x, :i ; r/2), j = (),...,/>} has a nonempty intersection. 

Cech complexes are higher-dimensional analogues of geometric graphs, a no¬ 
tion more familiar to probabilists. 

DEFINITION 1.2. Given a finite set X C W l and a real number r > 0, the 
geometric graph G(X. r ) is the undirected graph with vertex set X and edges 
[x, y] for all pairs x, y £ X for which \\x — y\\ < r. 

For a given a Cech complex, it is immediate from the definitions that its 1- 
skeleton is actually a geometric graph. Examples of both are given in Figure 1 . 



C(X.r) 


G(X.r) 


Fig 1. Take X = {an,..., xs} C R 2 . Note that the 2-simpIex [an, an, * 3 ] belongs to C(X, r ), since 
the three balls with radius r centred at X\,X 2 , X 3 have a common intersection. The Betti numbers of 
C(X,r) are /3 0 = 1 (one connected component), j3\ = 2 (two closed loops), while all others are 
zero. The geometric graph G(X, r) shows the 1-skeleton of C(X , r). 


A typical TDA paradigm is to create either U(X, r ) or C(X, r) as an estimate 
of some underlying sub-manifold, A4 C M. d , from which X is actually sampled, 
and then consider its homology, typically via what is known as persistent homology 
(with which we shall not concern ourselves in this paper) or through Betti numbers. 

We shall concern ourselves primarily with Betti numbers, a basic quantifier in 
Algebraic Topology. Roughly speaking, given a topological space X, the 0-th Betti 
number counts the number of its connected components, while for k > 
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1 , the k- th Betti number 8k(X) counts the number of A;-dimensional ‘holes’ or 
‘cycles’ in X. For example, a one-dimensional sphere - i.e. a circle - has 80 = 1, 
8 \ = 1, and 8 k = 0 for all k > 2. A two-dimensional sphere has 8 o = 1, 
8 \ = 0, and 82 = I - and all others zero. In the case of a two-dimensional torus, 
the non-zero Betti numbers are 80 = 1, 81 = 2, and 82 = 1- At a more formal 
level, 8 k(X) can be defined as the dimension of the k-tli homology group, however 
the essence of this paper can be captured without knowledge of homology theory. 
In the sequel, simply viewing 8 k(X) as the number of /.-dimensional holes will 
suffice. The readers who are interested in a thorough coverage of homology theory 
may refer to [15] or [27]. 

Returning to the two sets U(X,r/2) and C(X,r), by a classical result known 
as the Nerve Theorem ([9]) they are homotopy equivalent, and so from the point of 
view of TDA they are conceptually equivalent. However, since the definition of the 
Cech complex is essentially combinatorial, it is computationally more accessible, 
and so of more use in applications. Hence we shall concentrate on it from now on. 

There is now a substantial literature, with [18, 19] being the papers that moti¬ 
vated us, that shows that, given a nice enough M., and any <5 > 0, there are explicit 
conditions on n and r such that the homologies of U and C are equal to the homol¬ 
ogy of M with a probability of at least (1 — 5). Typically, these results hold when 
the sample X is either taken from M itself, or from A4 with small (e.g. Gaussian) 
random, perturbative error. However, if the error is allowed to become large (e.g. 
has a heavy tailed distribution) then these results fail, and a phenomenon known 
as crackle occurs. An example is given in Figure 2, taken from [1], in which the 
underlying Ad (a stratified manifold in this case) is an annulus in M 2 . From the 
point of view of TDA, the main aim of the current paper is to derive rigorous re¬ 
sults describing the distribution of the ‘extraneous homology elements’ appearing 
in Figure 2(d). 

These results, as central as they are to the paper, will be given in detail only 
in Section 5. There we shall show that the Betti numbers of the sets generating 
this spurious homology, which occurs far from the support of AT satisfy Poisson 
limit laws. The limit laws are stated in terms of the distance of these sets from 
Ad, so that they are in terms of limiting Poisson laws, the intensities of which 
depend significantly on the tail decay rate of the perturbative error. The practical 
importance of this is that it gives useful information on how much crackle occurs, 
what its topological nature is, and where it occurs, at least in a limiting scenario. 
The results of Section 5, therefore, are the main ones from the point of view of 
TDA. 

The reason that we need to wait so long to get to these results is that the natu¬ 
ral probabilistic setting for both stating and proving the theorems from which they 
follows is that of EVT. Recall that the classical setup of (multivariate) EVT takes a 
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Fig 2. (a) The original space M (an annulus) that we wish to recover from random samples, (b) 
With the appropriate choice of radius, we can easily recover the homology of the original space from 
the union of balls around a random sample from A4. (c) In the presence of bounded noise, homology 
recovery is undamaged, (d) In the presence of unbounded noise, many extraneous homology elements 
appear, and significantly interfere with homology recovery. 


point process generated by independently and identically distributed (i.i.d.) random 
variables in R^ - a random point cloud - scales and centers the cloud, and then, 
typically, shows weak convergence of the new cloud to a Poisson random mea¬ 
sure. The literature dealing with this setup is rich, and impossible to summarise 
here, but the main references for EVT are probably [14], [22], and [24], while a 
comprehensive introduction to point process theory is provided by [11], 

Beyond these more classical papers, over the last decade or so there have been 
a number of papers treating geometric descriptions of multivariate extremes from 
the perspective of point process theory, among them [4], [5], and [6]. In particular, 
Poisson limits of point processes possessing a U-statistic structure were investi¬ 
gated by [10] and [26], the latter also treating a number of examples in stochastic 
geometry. 

All of these papers focused on the (limiting) distributions of extreme sample 
clouds, but so far, surprisingly few attempts have been made at analyzing their 
topological features. The objective of the present work is to present a unified, ge¬ 
ometric, and topological approach to the point process convergence of extreme 
sample clouds. Results on the limiting behaviour of geometric complexes are quite 
new, with results not unrelated to our own, include [16], [17], [28], and [29]. In 
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particular [17] and [29] derived various central and Poisson limit theorems for the 
Betti numbers of the random Cech complexes C(X n . r n ), with X n a random point 
set in W 1 and r n a threshold radius. The resulting limit theorems depend heavily 
on the asymptotics of nrf, as n -A oc. For example, [17] investigated the sparse 
regime (i.e. nr^ —> 0) so that the spatial distribution of complexes is sparse, and 
they are observed mostly as isolated components. In contrast, the main focus of 
[29] was the thermodynamic regime (i.e. nr^ -A £ £ (0, oo)) in which complexes 
are large and highly connected. 

However, with the single exception of [1], already discussed above, none of 
these papers has results related to crackle. The main discovery of [1] was that the 
topological behaviour of the individual components of geometric complexes built 
over point clouds of (typically) heavy tailed random variables is determined by 
how far away they lie from A4. In the simple case in which Ad is a single point 
at the origin, this is exhibited by the formation of ‘rings’ around the origin, each 
ring containing points of the cloud which exhibit quite different geometric and 
topological properties. The positions of the rings depend crucially on the tail decay 
rate of the underlying probability distribution of the observations as well as the 
topological property being studied. The contribution of the present paper is to take 
the discovery of [1], made on the basis of only two examples, and develop it into a 
full theory, covering wide classes of distributions and developing far more detailed 
descriptions of the crackle phenomenon. 

The remainder of the paper is structured as follows: In Section 2, we establish 
general results on point process convergence under certain geometric constraints. 
Since the results here are very general, the conditions of the theorems may at first 
seem rather unnatural and hard to check. That this is actually not the case be¬ 
comes clear when we turn to specific scenarios depending on the tail properties of 
the underlying distributions. In doing so, it turns out to be natural to distinguish 
between what we shall henceforth call ‘heavy-tailed’ (e.g. power law tails) and 
‘light-tailed’ (e.g. e - ^^, r > 0) densities. Section 3 defines heavy-tailed densi¬ 
ties more precisely, via regular variation of tails (giving marginal distributions in 
the max-domain of attraction of the Frechet law) and then applies the general re¬ 
sults of Section 2 to obtain limiting Poisson results for this case. Section 4 then does 
the same for light-tailed densities (with tails controlled by the von Mises function). 
Typically, a standard argument in EVT characterises distributional tails via distri¬ 
bution functions (not densities), but our calculation will proceed relying heavily on 
the treatment of densities. Therefore, in both Sections 3 and 4, we shall distinguish 
the underlying distributions using density functions. 

The final Section 5 presents several applications of the point process conver¬ 
gence, one of which is the limit theory of Betti numbers alluded to above. In ad¬ 
dition, we give some independently interesting limit theorems on partial sums and 
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m a xima under geometric constraints. To the best of our knowledge, these results 
have no precursors in the random topology literature, since, while they are simple 
corollaries of the Poisson convergence that is our main theme, proving them ab 
initio would be quite hard. 

A final remark is that all random points in the present paper are generated from 
an inhomogeneous Poisson point process in W 1 with intensity nf (/ is a density 
function). However, all the results obtained can be carried over to the situation of an 
i.i.d. random sample without difficulty. Further, we only consider spherically sym¬ 
metric densities. Although the spherical symmetry assumption has very little to do 
with our results, we adopted it in order to avoid unnecessary technicalities, which 
would otherwise blur the message of the paper as well as making the treatment 
much longer and more cumbersome. 

Since the paper is rather long, and unavoidably heavy on notation, we conclude 
the Introduction with a short notation guide. 

2. Point Process Convergence: General Theory. In this section we prove 
a general convergence result for point processes with geometric constraints. The 
main result. Theorem 2.1, which holds under somewhat opaque conditions, will 
then be applied in the two following sections to obtain results under specific, and 
quite natural, distributional assumptions. Before getting to results, however, we 
need to set up a framework and notation, which we do in the following subsections. 


2.1. Point clouds and indicator functions for geometry. We start with a col¬ 
lection X = {Xi , i > 1} of i.i.d., Revalued random variables with spherically 
symmetric probability density /. For n > 1, let N n be a Poisson random variable 
with mean n, and independent of the X,. Then the random measure with atoms at 
the points X\,... , X^ n , which we denote by V n , is an inhomogeneous Poisson 
point process with intensity nf. 

Now choose a positive integer k, which will remain fixed for the remainder of 
the paper. Although some of our results hold for k = 1, we shall henceforth take 
k > 2, unless otherwise stated. As will soon be clear, in the k = 1 case many of the 
objects and functions we build will degenerate, in the sense that the corresponding 
results will become known results of non-topological EVT. 

Let 

(2.1) X k = {X u ...,X k } 


be the first k points and for an ordered fc-tuple i = (i\,..., if), 1 < i\ < ■ ■ ■ < i k , 
consider the /c-tuple of points 


= {X h ,...,X ik }. 


(2.2) 
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Symbol 

Definition or Meaning 

Place Defined 

C{X,r) 

Cech complex 

Definition 1.1 

G(X,r) 

Geometric graph 

Definition 1.2 

X 

Point cloud {A'i,.. . } 


Vn 

Poisson point process with intensity nf 


x k 

First k points in X 

(2.1) 

?C\ ,1 {i i 5 • • • j i k } 

Subset (Aij,..., Xi k } of X 

(2.2) 

2 - n,k 

All subsets of (1,..., n} of size k 

(2.3) 

Tm C-k,ni dk,ni Rk,n CtC. 

Scaling sequences, specific to section 



Xi, translated and scaled 

(2.10) 

4 n) 

First k points of normalised XiS 

(2.11) 

X i (n) ,i = {h,...,i k } 

Subset {X^\...,X^} 

(2.12) 

S{x) 

Maps ler onto +ve orthant of S d 1 

(2.13) 

L k 

Upper-diagonal subset of a fc-dim. cube 

(2.14) 

N^ k) 

No. of fc-tuples satisfying a condition 

(2.15) 

N^ k) 

No. of isolated fc-tuples as above 

(2.16) 

P(x;r) 

Local integral of density 

(3.7) 

m 

Polar lacobian 

(4.9) 

x A y, for x, y £ R d 

Xi <yi,i= 


x A y, for x, y £ R d 

Xi < yu i = 1 , ■ ■ ■, d 


a V b, for a, b £ 8 

max{a, b} 


a A b, for a, b £ R 

min{a, b} 


B(x;r) 

Ball, centre x, radius r 


INI 

Euclidean norm 


s d 

d-dimensional unit radius sphere 


Sd. 

Surface area of S d 


A k 

/c-dimensional Lebesgue measure 


Poi{\) 

Poisson distribution with mean A 


=> 

Weak convergence 


0 

Vector of zeros 


ei 

Vector with 1st comp. 1 and others zero 


Ann(a, b) 

Annulus, {a; £ R d : a < ||a:|| < b} 



We are interested in the geometric objects built over the Pi\, for all i in the set 
x \v n \,k> where 

(2.3) (J*i •>•’•') 'ik) • 1 — ^ ^ 'i'k — rnj-, 

and \V n \ = N n is the number of atoms of V n . 

To set this up formally, we define an indictor function h : (W l ) k —> {0,1} 
which will capture the geometry. Two concrete examples related to problems that 
will be considered in detail later sections are 

h{x i, ...,x k ) = l(G({xi,... ,x k }, 1) = T), 


(2.4) 
and 

(2.5) 


h(xi,...,x k ) = l(Pk -2 [C{{x t,...,x fc }, 1)) = l), 
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where 1 is the indicator function, G is a geometric graph, T is a fixed connected 
graph with k vertices, = denotes graph isomorphism, and C is a Cech complex. 
Throughout the paper we shall require that h is shift invariant, viz. 

(2.6) h(x i,...,Xk) = h(xi+y,...,xk + y), for all x±,..., x k , y e R d . 

We shall refer to this as the shift invariance condition in what follows. Further, we 
shall typically require a points in proximity condition, defined by the requirement 
that there exists a finite M > 0 for which 

(2.7) h(0,x\, ...,Xfc-i) = 0 if ||xj|| > M for some i E {1,... ,k — 1}. 

In addition, given anon-increasing sequence of positive numbers (r n , n > 1) (the 
case r n = constant is permissible) we define scaled versions of h, h n : (M.' r ) k —>• 
{0,1}, by setting 

(2.8) h n (xi,...,x k ) = h(xi/r n ,... ,x k /r n ). 

Clearly, h n is also shift invariant. Furthermore, if r n is small, then (2.7) implies that 
h n (x i,..., Xk) = 1 only when all the points x \,..., Xk are close to one another. 

In what follows, we shall be particularly interested in A:-tuplcs X\, i G V n \,k 
which not only satisfy the constraint implicit in the indicator h n but are also sep¬ 
arated (by at least r n ) from the other points in V n . Thus, we need to introduce 
another sequence of indicator functions. Given a subset X of k points in W 1 , and a 
larger, but finite, set 3^ D X of points in the same space, define 

(2.9) 

g n (X, y) = h n (X) x 1 (G(X, r n ) is an isolated component of G(y, r n )). 

2.2. Sequences of point processes. Since we are going to be dealing with the 
geometry of random variables ‘in the tail’ of /, as in standard EVT we are going 
to need to normalise them before trying to formulate any results. 

Thus, with k still fixed, we assume that we have two sequences of constants, 
dk,n £ R and Ck. n > 0, with which to shift and then scale the X,, and, in notations 
that will save a lot of space, define: 

(2 IQ) ^(n) A Xj — dk,nS(Xj) 

1 C k ,n 

(2.11) X^ n) = {x{ n) ,...,xf} 

(2.12) X[ n) = {X^,...,X^} 
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where 

(2.13) 


S(x) = 



maps points x = (.x 1 - 1 ),..., x (,l> ) in W l onto the positive orthant of S d_1 . 

We can now define the two random measures, that will be at the centre of all that 
follows, over the space E' k x L k , where E k is a measurable set in ([—oo, oo] d ) 
and L k is the upper-diagonal subset of a /.--dimensional cube, viz. 


(2.14) L k = {(z 1 ,...,z k )e[0,l} k :0<z 1 <---<z k <l}. 


The measures, for Borel A C E k x L k , and e x a Dirac measure at x £ E k x L k , 
are given by 

(2.15) NW(A) i Y h -( X 0e (x M ./\ Pn]) (A), 

1&X \ Vn\,k 

and 

(2.16) NW(A) ± Y 9n(Xi,V n )e {x ( n) . /]rnl) (A). 

' 1 &X \V n \,k 

(Note that = 0 whenever k > n.) 

(k) 

The point process Nn ; counts A;-tuplcs satisfying the geometric condition im- 

~(]A 

plicit in h n , while A r „ ' adds the additional restriction that the geometric graphs 
generated by such /c-tuples be distance at least r n from all the remaining points 
of V n . It will turn out that, while A„ always dominates N' n , both exhibit very 
similar limiting behaviour. In fact, the proof of the general result on point process 
convergence will follow the route of first establishing a result for Nn , which is 
somewhat easier because the geometric condition implicit in the indicator does not 
include V n , and then showing that since the two processes are close to one another 
the same limit holds for Nn . 


2.3. A general theorem. Retaining the notation of the previous subsection, we 

now set up the main conditions that will ensure the weak convergence of the pro- 

(k) ~(k) ^ 

cesses Nn and Nn . Although we ultimately want theorems in which conditions 

relate to the tail behaviour of the underlying density of the X. t , the conditions in this 
section relate both to this and to what amounts to asymptotic sparsity for certain 
structures in V n . In Sections 3 and 4 we shall see how to derive this independence 
from tail considerations. 
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For our first condition, which we refer to as the vague convergence condi¬ 
tion, we assume that there exist a Radon measure u k on a measurable set E k c 
([—oo, oc] d ) k , not assigning mass to points at infinity, and normalising constants, 
as in the previous subsection, such that 

(2.17) (fj p [h n {x k ) = l, 4 n) e •} A !/*(■), 

where A denotes vague convergence in E k . 

The second, sparsity condition, is given by 

(2.18) 

n 2k ~ l P {h n (X k ) = h n (X !,..., X t ,X k+1 ,..., x 2k _ l ) = 1, X^_, € K } -> 0 , 

as n —> oo, for every l = 1,..., k — 1 and every compact set K in E 2k -i■ This 
condition prevents the random /.--tuples from clustering too much, and will lead to 
a Poisson limit for iVA . The same kind of sparsity, or anti-clustering, condition 
can be found in, for example, [10] and [26]. 

~(k) 

An additional condition will be needed for the convergence of NX . For this, 
let V' n be an independent copy of V n . Then this condition, combined with (2.17), 
requires that, asymptotically in n, the k points in X k all be distance at least r n from 
the points in V' n . Specifically, 

(2.19) Qp {g n (X k , X k UV' n ) = l, xf, n) € •} 4 */*(■)■ 

As one might expect from these conditions, when they hold the weak limits of 

(k) ~ (k) 

Nn and A r A will typically be expressible via Poisson random measures, which 
we denote by N- k \ and now define. 

Writing X k for /.--dimensional Lebesgue measure, the Poisson random measure 
N- k> on E k x with mean measure u k x \ k (with X k restricted to L k ) is defined 
by the finite dimensional distributions 

P (N^ k \A) =m ) = e ~^ kXXk ^ A \(iy k x X k ){A)) m /m\, m = 0,1,2,..., 

for all measurable A C E k x L k with (v k x X k )(A) < oo. Furthermore, if 
Ai,... ,A m are disjoint subsets in E k x L k , then N^ k \A\),... ,N^ k \A m ) are 
independent. 

The random measures N^ k \ Nn k \ and A r - fc / are all regarded as random elements 
in the space M p (E k x L k ) of Radon point measures in E k x L k . Details on such 
random measures, including issues related to their weak convergence in the space 
of Radon measures, along with the vague convergence of (2.17) and (2.19), can be 
found, for example, in [22] and [23]. 
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Before stating the main result of this section, a few words about the above three 
conditions are in order. Firstly, we note that more specific definitions of the state 
space Ek and concrete definitions of the normalising constants (c^ n , dk. n ) will 
be given in the subsequent sections, where they will be seen to be dependent on 
whether the underlying density / has heavy tail or light tail. 

In the heavy-tailed density case, we take Ek = ([—oo, oo] d ) fc \ {0} and c4 n = 
0. Thus there is no shift of the points Xi by ,S’(X,), and all the scaling does is 
bring points closer to the origin. Consequently the point processes Nn and Nn 
asymptotically count geometric events equally likely in all the orthants. 

On the other hand, in the case of light-tailed densities, one needs to take Ek = 
((—oo, oo] d ) k and c4 n is generally non-zero. Since X r is translated by dk %n S(Xi), 
this implies that, at least when n is large, Nn and Nn count geometric events 
basically only when they occur in the non-negative orthant. Nevertheless, the as¬ 
sumed spherical symmetry of the density allows our main results to be extended to 
general domains with only minor additional arguments. 

Thus, since the aim of this section is to establish general theory on point process 
convergence, we shall leave this interperative complication until it arises later in 
specific examples. 

As remarked above, throughout the paper we take k > 2. Given the definitions 
above, we can now explain why this is the case. If k = 1, there is no shift invariant 
indicator to associate with the point process, and we have 



This type of the point process has been well studied in classical EVT, and we have 
nothing to add about it. 

We can now state the main result of this section. The proof is deferred to the 
Appendix. 


THEOREM 2.1. Let (Xi, i > 1 ) be a sequence ofi.i.d. W 1 -valued spherically 
symmetric random variables. Let h n : (W l ) k —> {0,1} be a sequence of indica¬ 
tors as in (2.8), satisfying shift invariance as in (2.6) and the points in proximity 
condition (2.7). Assume also that the vague convergence condition (2.17) and the 
sparcity condition (2.18) hold. Then 

N (k) ^ jV (fc) j 

where =e denotes weak convergence in the space M p (Ek x Lk). 

In addition, suppose that the global separation condition (2.19) holds. Then we 
also have that 


=> N {k) . 
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REMARK 2.2. Note that the case h = 1 is precluded by (2.7). In fact, ifh = 1, 
then it is easy to see that (2.17) and (2.18) cannot hold simultaneously. 


3. Heavy Tail Cases. As a special case of the general theory introduced in 
the last section, we shall discuss point process convergence when the underlying 
density on U. d has a heavy tail. More specifically, we assume that the density has 
a regularly varying tail (at infinity) in the sense that, for every 6 £ S' 1 ' 1 (equiva¬ 
lently, for some 6 £ S d ~ l , due to the spherical symmetry of /), and some a > d, 

f (tvO^) 

lim = t ~ a , for every f > 0. 

r-roo j{r0) 

Writing RV- a to denote the family of regularly varying functions (at infinity) of 
exponent —a, then the above is equivalent to the requirement 

(3.1) / G RV- a . 


In the one-dimensional case (d = 1) it is known that regular variation of the tail in 
/ suffices for the distribution to be in the max-domain of attraction of the Frechet 
law; see for example Theorem 3.3.7 in [14]. 

In general, given a non-increasing sequence of positive numbers (r n , n > 1), 
we define h n , g n as in Section 2. Then, for a fixed positive integer k > 2, the point 
processes we are going to explore are 


(3.2) 


and 

(3.3) 


Ni k \.) 


- hn ^ e (R^, i/|Tn|)(')’ 

l£X \V n \,k 



^ Sn(*i,P„)e( ^ i/| Pn |)(-) 

lGX \Vn\,k 


where , n > 1) is asymptotically determined by 

(3.4) n k r d ^RlJ(R ktnei ) k ^ 1 as n ^ oo , 

with e\ = (1,0,..., 0)'. 


Remark 3.1 . Note that (3.4) implicitly precludes very fast decay ofr n to zero. 
In fact, if n r, i -A 0 as n —» oo, then (3.4) implies that we must also have 
that R d n f{Rk,n&i) k —>• oo, which contradicts (3.1). However, under the implicit 
restriction that such rapid convergence does not occur, (3.4) effectively determines 
the Rk, n as an intrinsic function ofk and r n . 

Note also that the case r n = 1 was the one treated in [1], although for a far 
more limited class of densities. In the case r n = 1, it is easy to check that Rk, n G 
RV^u-d/k)- 1 as a function of n. 
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The convergence in law of and takes place in the space M p ( E k x 
L k ), where E k = ([— oo, oo] rf ) fc \ {0} and L k is given in (2.14). Let N^ be a 
limiting Poisson random measure, whose intensity is, as in Section 2, denoted by 
i-'fc x Afc, where Xk is the A;-dimensional Lebesgue measure concentrated on L k . As 
a consequence of the heavy tail of /, the measure v k also exhibits a power-law tail 
structure. More specifically, for a rectangle (ao,&o] x ••• x (a /,--1 ■ b k — i1 C Ek = 
([-00,00] d )*' \ {0}, which is bounded away from the origin, 

(3.5) ^fc((a 0 ,& 0 ] x ••• x (a k -i,b k -i\) 

= 77 / h(0,y)dy [ \\x\\- ak dx 

,/(R d ) fc—1 Jai^x^bi , i= 0 ,...,k —1 

where -< and A indicate componentwise inequalities. 

It is worth mentioning that v k exhibits the scaling property: 

(3.6) u k {sA) = s-^ k ~ d ^ k (A) 

for all s > 0 and measurable A. We shall also need one additional function. For 

x = {x\,... , x k ) €E (M. d ) k , r > 0, we define 

(3.7) p(x;r) = [ f(z)dz. 

We now have all we need to formulate the following result. The proof is given 
in the Appendix. 

THEOREM 3.2. Let (Xj, i > 1 ) be a sequence ofi.i.d. W 1 -valued spherically 
symmetric random variables with density f having regularly varying tail as in 
(3.1). Let h n : (R d ) fc —> {0,1} be a sequence of indicators as in (2.8), satisfying 
the shift invariance of (2.6) and the points in proximity condition (2.7). Then, the 
point processes N^ and given by (3.2) and (3.3) weakly converge to N^ 
in the space M p (E k x L k ). 

To make the implications of the theorem a little more transparent, consider the 
simple case for which 

( 3- 8) f(x) = ° „ Q , 

1 + ||x||“ 

and assume that r n = n s with —k/d(k — 1) < s < 0. Then, (3.4) reduces to 
n k+sd(k-i) (jk ftd-ak ^ an d so i v j n g this with respect to R k , n gives 

R k n = (ji^-d/k)- 1 ^[l+sdil-k^j/la-d/k) ^ 
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0 


\ 

) 


#r 2 ' 


\ #r 2 ■ 


#r fc -x ~ oo #r fc _ 2 ~ oo 
#r». ~ Poisson / #r fc _! ~ Poisson 
#r fc+1 ~ 0 / #T k ~ 0 



Poisson 

0 

0 


Fig 3. Layer structure of annuli for random geometric graphs. For the density (3.8), Rk,n is a 
regularly varying sequence with exponent [1 + sd(l — fc -1 )]/(a — d/k). The number of subgraphs 
isomorphic toYk outside B(0; Rk.n ) is approximately Poisson. 


which in turn implies • • • << R kn « R k -i, n <<•••<< R2,n- 

Being more specific on the geometric side, consider the geometric graph exam¬ 
ple of (2.4), where connected graphs Twith k vertices are fixed for k = 2,3,.... 
Then the theorem implies that W 1 can be divided into annuli of increasing radii, 
and the Tfc, A: = 2,3,... are (asymptotically) distributed in these annuli in a very 
specific fashion. Letting Ann (A. L) denote an annulus of outer radius L and inner 
radius K, we have (in an asymptotic sense) 

• Inside Ann(f? 2 ,n, oo), there are finitely many I^, but none of T 3 , T^.... 

• Inside Ann(f? 3 jn , f? 2 , n ), there are infinitely many T 2 and finitely many r^, 
but none of T^ Ts,.... 

In general, 

• Inside Ann(R kjU , R k ~ i, n ), there are infinitely many T 2 ,..., Tfc_i and finitely 
many T k , but none of r fc+ 1 ,r fe+2 , .... 


4. Light Tail Cases. This section treats point process convergence when the 
underlying density on W 1 possesses a relatively lighter tail than in the previous 
section. Typically, in the spirit of extreme value theory, light-tailed densities can be 
formulated by the so-called von Mises function. In particular, in a one-dimensional 
case (d = 1), the von Mises function plays a decisive role in a characterisation of 
the max-domain of attraction of the Gumbel law. See Proposition 1.4 in [22]. The 
density formulation similar to our setup can be found, for example in [3] and [4]. 

Our results for the light tailed case fall into two categories. In the first, we have 
results more or less paralleling those of the previous section, although the normali¬ 
sations and limiting distributions are somewhat different. Nevertheless, the ‘annuli 
structure’ in the preceding section carries through. This case is treated in Section 
4.1. Recall, however, that in the heavy tail case these annuli are rather thick, espe¬ 
cially when compared the typical size of the geometric objects, which is of order 






POINT PROCESSES UNDER GEOMETRIC CONSTRAINTS 


15 


r n . In the light tail scenario, this is not the case, and the annuli are actually quite 
thin, which leads to questions that are not present in the heavy tail case. These 
questions are posed and solved in Section 4.2. 

4.1. Annuli. The light tail case considers (Xj, i > 1) of R d - valued i.i.d. spher¬ 
ically symmetric random variables with density given by 

(4.1) f{x) = L(||x||) exp{-^(||x||)} , x€R d . 

Here, ip : M + —> R is a C 2 function of von Mises type, in that 

ip'(z) > 0, ip(z) — > oo, (l /ip'}' {z) 0 


as z —> Zoo, for some € (0, oo]. The main objects of this paper are the geo¬ 
metric objects far away from the origin, so the bounded support of the density (i.e., 
Zoo < oo) is relatively less interesting to us. Therefore, in what follows, we restrict 
to the unbounded support case, in which = oo. For notational convenience, we 
introduce the function 


Since a!{z) —> 0 as z —> oo, the Cesaro mean of a' converges as well. That is, 

Z 

a'(r)dr -A 0 , as z —> oo . 

Suppose that the function L : l i -a R f in (4.1) is flat for a; viz. 


ajz) =- ( 

z z Jo 


(4.2) 


L(t + a(t)v) 

-—--—> I as t —> oo uniformly on bounded v-sets. 

Lit ) 


This assumption implies that L is negligible in its tail, and hence the tail behaviour 
of / is determined only by ip. Here, we need to introduce an extra regularity con¬ 
dition on L. Assume that there exist 7 > 0, zq > 0, and C > 1 such that 


(4.3) 


L(zt ) 

W) 


<Cf for all t > 1 , 2 > zq . 


In view of Corollary 2.0.5 in [7] (or Theorem 2.0.1 there), 


(4.4) 


lim sup 

Z—>-00 


L(zt) 

L{z) 


< 00 for all t > 1 


suffices for (4.3). Observe that if L is a polynomial function, then both (4.2) and 
(4.4) are satisfied. On the other hand, (4.4) does not hold if L grows exponentially. 
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Given a non-increasing sequence (r n , n > 1) of positive numbers (as usual con¬ 
stant r n is permissible), we define h n and g n as in the previous sections. Suppose 
that there exists a non-decreasing sequence (R k , n , n > 1) determined by 

(4.5) n k r^ k ~ l) a(R k , n )R^ f(R k>n e i) fc -5- 1, n -A oo . 

As pointed out in Remark 3.1, note that arapid decay of (r n ), e.g., n k r^ k 1 ' ) —> 0, 
n -A oo is not permissible. For k > 2, using abbreviations (2.10), (2.11), and 
(2.12), so that 

Cfc,n Rk,n ) j dk,n Rk,n i 

we shall consider the following point processes in the space M p (E k x Ly.) with 
E k = ((— oo,oo] d ) k and L k given by (2.14): 



and 


(4.7) 



Y 9 n(X { ,V n ) 

i£l \Vn\,k 


e (x.^\i/\r n \)^ 


(k) ~(k) 

As seen below in more detail, the nontriviality of weak limits of Nn and Nn 

mainly depends on the limit value of a(Rk : n)/^n- Indeed, if a(Rk, n )/fn h as a 
(k) ~(k) ’ 

non-zero limit, then iV„ and Nn both converge to the same Poisson random 

measure. On the contrary, if a{R.y. n ) An -A 0, then Nn } goes to 0 in probability 

~(k) 

(and of course, so does Nn ), which implies that the geometric objects implicit 
in h n cannot be observed outside 77(0; Rk,n), regardless of whether or not those 
objects are isolated from other points. 

The following is the main limit theorem of this subsection. The proof is deferred 
to the Appendix. 


THEOREM 4.1. Let (2Q, i > 1 ) be a sequence ofi.i.d. W l -valued spherically 
symmetric random variables with density f provided in (4.1). Let h n : (M rf ) fc —» 
{0,1} be a sequence of indicators as in (2.8), satisfying the shift invariance of (2.6) 
and the points in proximity condition (2.7). 

( i ) If a(Rk : n)/ r n -A-cG (0, oo] as n —> oo, then the point processes N^ and 
Nn given by (4.6) and (4.7) converge weakly to a Poisson random measure with 
intensity u k x X k . Here, Ay. is the k-dimensional Lebesgue measure concentrated 
on L k , and the measure u k is given by 


(4.8) v k ((a 0 ,b 0 ] x ••• x (a k -i,b k -i\) 
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x 1 (a,i A (p + c~ 1 (9,y i ))e i = l,...,k — lj 
k -1 

x exp{ -/qo - cT 1 ^(<9, z/j)} /z(0, y) dpJ(6)dddy, 
i= 1 

where <n, bi, i = 0..... A: — 1 are d-dimensional real vectors with — oo -< a, A 
bi A oo, and J(9) = \dx/dd\ is the Jacobian 

(4.9) J(0) = sin fe - 2 (0 1 ) sin fc_3 (0 2 ) ... sin(0 fe _ 2 ). 

(zz) Suppose a,{Rk,n)/ r n —>• 0 as n —> oo a/z<7 a zs eventually non-increasing (i.e., 
a is non-increasing on (.xq, oo)for some large xq). Then, the point processes Nn ' ! 
and Nn ' converge to 0 in probability. 


As was the case in the heavy tailed scenario, Theorem 4.1 (i) again implies the 
layer structure of annuli at different radii • • • << 74, n « 74_i,n <<•••« 
7?2,n in which different structures can be found. In the current case, however, this 
is not as immediate, since now the normalisation also involves a translation, the 
asymptotic effect of which is to count geometric events only when they occur in 
the nonnegative orthant. 

Recall that we are particularly interested in the distribution of geometric events 
outside the ball 77(0; Rk, n )i k = 2,3,..., and this is captured by 

^2 9n(X\,V n ) 1(1174. || > 74,n, j = l,...,k) . 

' &X \V n \,k 


However, by the spherical symmetry of the underlying density, along with the in¬ 
dependence of the sample points, the weak limit of this quantity is the same as that 
of 

2 d Y, 9n(X i,Pn)e Ar Cn)(([0,oo] (i ) fc ), 

ieX \v n \,k 

the limit of which can be computed directly from Theorem 4.1 (?'). In other words, 
the spherical symmetry of a density allows one to extend the basic result on the 
weak convergence occurring in the nonnegative orthant to that in all orthants. Con¬ 
sequently we can obtain the same qualitative separation of geometric events into 
distinct annuli that we saw in the heavy tailed case. 

To see how this works, consider the simple example f(x) = C'e H x l7/ r , for 
some 0 < t < 1, and take r n = 1. Then, a(z) = z 1 ~ T clearly has a non-zero limit 
as 2 —> 00 , and so Theorem 4.1 (z) applies. Then, (4.5) becomes 


n 


k jjd-r (jk e ~ kR l. Jt 


18 


OWADA AND ADLER 


and the solution 

(4.10) Rk,n = (t log n + AT 1 (d - r) log(r log n) + r log C) 

grows only logarithmically in n, whereas the /?/,. „ of the previous section exhib¬ 
ited, essentially, power law growth. Thus, the description inherent in Figure 3 re¬ 
mains unchanged, except for the change in the values of Rk, n - 

REMARK 4.2. In particular, if a(Rk n )/r n —> oo, the limiting intensity mea¬ 
sure Vk can be simplified to 

(4.11) u k ((a 0 ,b 0 \ x ••• x (a fc _i,6 fe _i]) 

= 77 f h(0,y)dy I f e~ kp dpj(9)d9 , 

K - J( R d ) fc - 1 J Jai^pO^bi, *=0,...,fc— 1, 0X0 

where a,, b r , i =- 0..... A' — 1 are d-dimensional vectors such that — oo -< a* ^ 
bi R oo. 

REMARK 4.3. 7/2 f/ie cave ofk= 1, //ie point process Nn ! is no longer asso¬ 

ciated with the indicator h n and is given by 



Interestingly, a simpler argument than the proof of Theorem 4.1 (or combining 
standard arguments of point process theory - e.g. Chapters 5 and 6 of [23] - with 
the Palm theory in the Appendix) shows that regardless of whether (i(R\, n ) has a 
zero or a non-zero limit, Nn 1 converges weakly to a Poisson random measure with 
a non-trivial intensity. 

Before concluding this subsection, we look at several examples of light-tailed 
densities for which the corresponding point processes exhibit different limiting 
behaviours. 

EXAMPLE 4.4. Suppose, for simplicity, that r n = 1 and consider the following 
probability densities on M d : 

fi(x) = Li(||x||)e" IN|T , 0 < r < 1, x€R d , 
h{x) = L 2 (||x||)e- ||a:||loglog||a: l |/log||x ll, x e R d , 

where the L\ and L 2 satisfy (4.2) and (4.3). The densities fi and / 2 are usually 
referred to as subexponential densities, since their tails decay more slowly than 
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that of an exponential distribution. For both f\ and fi, it is easy to check that 
a(z ) —> oo as z —> oo, so, by Theorem 4.1 (i), the point processes and 

weakly converge to a Poisson random measure with intensity ip. x A/,., where ip x - is 
given by (4.11). On the other hand, if the density has the same tail as an exponential 
distribution; for example, 

f 3 (x) = L 3 (\\x\\)e~ M , x£R d , 

then ip(z) = z and a(z) = 1. In this case, Nand once again weakly 
converge to a Poisson random measure. However, its intensity measure u;- x A/,, is 
more complicated, where ip, : is given by (4.8) with c = 1. We shall also consider the 
densities with more rapidly decaying tails than an exponential distribution (they are 
sometimes called superexponential densities). Two examples of superexponential 
densities are 


h(x) = L 4 (||x||)e l|:r||T , r > 1, x€R d , 
f 5 (x) = L 5 (||x||)e- ||a:||log||x|l/loglog||x|1 , xeR d . 

For fi and f§, it follows that a(z ) — > 0 as z — > oo, and so Theorem 4.1 (ii) implies 
that the point process Nf ' goes to 0 in probability. The densities fi, fi, and fr, 
differ only slightly in their tail behaviours, but the corresponding point processes 
possess totally different limits. Finally, we point out that even for fi and fr„ if one 
chooses (r n , n > 1) so that a{Rk,n)/ r n —> c £ (0, oo], then the point process can 
converge to a non-trivial Poisson random measure. 

4.2. At the annuli boundaries. The claim of Part (ii) of Theorem 4.1 is that, 
when a(Rk, n )/fn —■► 0, the geometric objects being counted do not exist outside 
of the ball i?(0; Rk /n f at l easl from the view of the point process convergence. We 
next want to explore the existence of the same objects inside B( 0; Rk,n), under the 
condition that these objects must be isolated from other random points by at least 
r n . This question was partially and negatively answered in [1], in the framework of 
the asymptotics of the expected Betti numbers of the Cech complexes associated 
with a random sample and a unit radius. In particular it was shown there that for 
the standard Gaussian distribution, all the expected Betti numbers of order k > 1 
vanish and the resulting Cech complex becomes contractible. In what follows, we 
continue working on the same question under the conditions of Theorem 4.1 (ii) 
from a more comprehensive viewpoint. 

Theorem 4.5 establishes the existence of two sequences of balls with different 
radii, the smaller ones ultimately containing so many points that they can be cov¬ 
ered by a union of balls with radius r n and centred on the points. On the other 
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hand, ultimately there are no points outside a larger balls. The main point, how¬ 
ever, is that the differences between the radii of the two balls decays at the rate of 
o(r n ). We conclude, therefore, that, at least when n is large, a union of balls with 
radius r n centred at points in V n becomes contractible, and accordingly, the geo¬ 
metric objects that we have been studying up until now, which are isolated from 
other points, fail to exist anywhere in all of R d . 

In order to get a clear picture, we shall add extra assumptions to the setup of 
Theorem 4.1 (ii). In particular, we assume that 


(4.12) 

ip £ RV v for some 1 < v < oo , 

(4.13) 

L = C (suitable normalising constant) 

and 



(4.14) (r n ) is a regularly varying sequence which decreases to 0 as n —>■ oo . 

The reason why we need v > 1 in (4.12) is as follows. By Proposition 2.5 in 
[23], which establishes the regular variation of the derivative of a regularly varying 
function, we now have 


a(z) = 1 /ip'(z) € RV i_„ . 

The setup of Theorem 4.1 (ii) requires that a(z) — > 0 as z — > oo, so the regular 
variation exponent of ip cannot be less than 1. Finally we observe that ip*~ e RV\/ v 
and a o ijX £ RVy_ v y v , and both functions are eventually monotone. 

THEOREM 4.5. Assume the conditions of Theorem 4.1 (ii), as well as (4.12), 
(4.13), and (4.14). Furthermore, assume that 

a o ip < ~( log n) 

(4.15) -- log log n -A 0 as n —>• oo . 

r n 

Then, there exist two sequences (Rn \ n > 1) and (Rn \ n > 1) 
n —> oo, 

p| B(0 ; rM) C [j B(X\r n ) , V n n B(0; R^)° = 0 
^ xeP„nH(o ; i4 0) ) 


such that, as 



and 


(4.16) 


^(r^-rW) ^ 0 . 
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Fig 4. Asymptotically, the smaller ball of radius is covered by a union of balls with radius 
r n centred around the points in V n . There are no points outside a larger ball with radius Rf\ The 
difference between Rn' 1 and Rn^ vanishes at the rate of o{r n ). 


REMARK 4.6. Typically, the solution to (4.5) is given by Rk, n = (log n + 
bk , n )» where bk, n / log n —» 0 as n -A oo; (cf. (4.10).) Because of the regular vari¬ 
ation of a o the condition a(Rk, n ) / r n -A 0 in Theorem 4.1 (ii) is equivalent 
to a o i/;' s_ (log n)/r n -A 0. In essence, therefore, assuming (4.15) adds a stronger 
condition to Theorem 4.1 (ii). We have not been able to determine whether or not 
the same result can be obtained when a o (log n)/r n -A 0 holds but (4.15) is 
no longer true. 

PROOF. First of all, we claim that (r n ) G RVq. To see this, note that it will 
suffice to show that log r n / log n —> 0 as n —> oo, since (r n ) is a regularly varying 
sequence, (cf. Proposition 2.6 (i) in [23].) Since a o G RV^_ v y v , it is clear 
that a o (log n) log n -A oo, and so, 

log r n a o (log n) r n \ogr n 

1 — TZTTj Tj ^ 0, n } oo. 

iogn r n a o ^(logn) logn 

For positive numbers g and p, let Q n {g, p) be a family of cubes with grid gr n 
that are contained in /i(0; p). Fix g > 0 sufficiently small so that 

{ Q n Pn 7 ^ 0 for all Q G Q n (g, p) } C | B(0; p) C \J B(X;r n )} 

X (z'Pn^B (0] p) 

for all p > 0 and n > 1. We define R'n ! and Rn' 1 as follows: 

R^ ] = (A n ), A n = log n + d log r n - log log r ~(log n) - 5, 


where <5 is a positive constant such that 
(4.17) d — e s g d C < 0 , 
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and 

RW=^(B n ), 

where 

B n = log n + (d — 1) log V’'* - (log n) + log a o ij)*~ (log n ) + log log n. 
We need to prove that, as n —> oo, 

(4.18) P{ Q n V n = 0 for some Q € Q„(g, ^r, 0) ) } ->■ 0 , 
and 

(4.19) p|p n n^(O;«i 1) ) c = 0} 1. 

The probability in (4.18) is estimated from above by 

y t P{0 n P n = 0} = ^ exp|-n j f(x)dx j 

QeQ»(s,^ 0) ) Q&Qn(g,Rn^) Q 

(4.20) < J] exp{-n(5r n ) d /(^ 0) ei)} 

QeQ„(g,fl£ 0) ) 

d 

exp{-g d nr^f(R^ei)}. 

By virtue of the inequality < i/W (log n) and (4.17), we have, as n —> oo, 

d\ogr~ l R^ - g d nr d f(R ( £ ) e i) < (d - e 5 g d C) logr“V^(logn) -> -oo 

from which the rightmost term in (4.20) vanishes as n —> oo. 

Next, we turn to proving (4.19). Since 

p{ Vn n B( 0; R^y = 0 } = exp j -n j (i) f(x)dx \ , 


it is enough to show n f, .. p) f(x)dx —> 0, as n —> 00 . By the polar coordinate 

0 |F|| — -^n 

transform with J(9) = |<9a;/<90|, we can write 
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where Sd-i is a surface area of the (d — 1)-dimensional unit sphere in W 1 . The 
dominated convergence theorem guarantees that the integral above converges to 
/o°° e ~ p dp = 1- Therefore, we only have to verify that 

na(R^)(R^) d 1 e~d( Rn ) —>• 0 , as n —>■ oo . 

Substituting Rn = ip*~(B n ), we have 




aoip^(B n ) ( ip^(B n ) 
a o ip*-(log n) \^d~(logn) 


d -1 

(log n) -1 . 


Since B n / log n —> 1, it follows from the uniform convergence of regularly varying 
functions (cf. Proposition 2.4 in [23]) that 

aoip^(B n ) ; 1 ip^(B n ) ; 1 
a o -ip^ (log n) ’ ip*— (log n) 


So the proof of (4.19) is complete. 

It remains to establish (4.16). The mean value theorem yields 

- ^i 0) ) = r~ 1 ('ip^)'(t n ) (B n - A„,) 

aoip^(t n ) a o ip*~ (log n) 

— ii-n \ Air), 

a o ip^ (log n) r n 

where t n lies in between A n and B n . Since A n / log n —> 1 and B n / log n — > 1, 
we have t n j log n —> 1, and thus, 


a o ip^(t n ) 


1 as n —> oo . 


a o 'i/d - (log n) 

To finish the argument, we have to establish the following three limits: 


(4.21) 

a o ip^ (log n ) 

log r n -A' 0, 

r n 

(4.22) 

a o ip^(log n ) 

log a o ip^ (log n) —> 0 , 

r n 

(4.23) 

a o ip^(log n ) 

log ip^ (log n) —>■ 0 . 

r n 

Since logo 

o ip*^ (log n) < log r n 

for sufficiently large n. 


(4.22). By virtue of the regular variation of a o r ijB and ip^, we have that 


log a o tp < ~(log n) 
log log n 


and 


/ log 'i/d - (log n) 
\ log log n 


are bounded sequences. Both (4.22) and (4.23) now follow from (4.15), and so we 
are done. □ 
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5. Applications. 

5.1. Limit Theorems for Betti Numbers. The results of the previous two sec¬ 
tions show the existence of a sequences of annuli containing different kinds of 
geometric objects. This subsection will further examine this layer structure and the 
topological properties of the objects they include, relying on the notion of Betti 
numbers to quantify the topology. Our aim is to derive limit theorems for the Betti 
numbers of the Cech complex built over V n = {X \...., Xy n }. where ( X r . i > 1) 
is an i.i.d sample drawn from a spherically symmetric distribution, and N n = \V n \ 
is a Poisson random variable with mean n and is independent of (X t . i > 1). For 
k > 3, we take h n : (M. d ) k -A {0,1} as in (2.5); viz. 

h n {x i, ...,x k ) = l(/3 fc _ 2 ... ,x k },r n )) = 1) , xi,...,x k e R d , 

and define g n : ( R d ) k -> {0,1} as in (2.9). Also, we define 

Sk,n= 9n(Xi,'Pn)l{\\X ij \\> R k , n , j = k) 

iGl \r„\,k 

for Aj = (X n ,...,X ik ) with i = (n,...,i k ) £ T\v n \, k - The definition of R, k>n 
depends on whether the underlying density has a heavy tail or a light tail. Specif¬ 
ically, if the underlying density has a regularly varying tail as in (3.1), then (3.4) 
defines the R k}n , while (4.5) determines the R k , n if the density is given by (4.1). 
We are interested in the behaviour of the Betti numbers 

P k - 2 (c{V n n B(0\R ktn ) c ,r n )y 

The most relevant study to this subsection is [1], in which the asymptotics of 
the expected Betti numbers were discussed. We, however, go well beyond this, by 
establishing Poisson limits for the Betti numbers. As originally given in [1], we 
shall provide one useful inequality to elucidate the relation between S k)U above 
and the Betti numbers: 

S k ,n ^ fik- 2 (c(V n n -B(0; R k ,n) C i r n)^ < S k ^ n + L k ^ n , 

where 

L k ,n= ^2 h n {Xi)l(\\X ij \\>R kj n,j = l,---,k + l), 

{&x \v n \,k+l 

with 

hn{x 1 ,.. ■ ,x k+1 ) = l(C({xi,.. .,x k+1 },r n ) is connected). 

The limit theorem below demonstrates that L k:n tends to zero in probability, and 
as a result, S k)U and /3fc_2 {C[V n FI B{ 0; R k , n ) c , r n )) asymptotically coincide. 
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Fig 5. Fork = 3, d = 2, the Betti number /3iyC (V„nB(0; R, 3 ,n) c ,r n )j counts (one-dimensional) 
cycles outside of B( 0; Ri, n ), while ignoring cycles inside the ball (e.g. (a), (b), and (c)). 


THEOREM 5.1. Under the assumptions and notations of Theorem 3.2, 
h- 2 [c{V n UB{Q-R Kn )fr n )) => Poi ^°’ y ) dy 

where Sd -1 is the surface area of the (d — l)-dimensional unit sphere in Ml 1 . 

PROOF. First of all, note that the variable Sf.. n defined above can also be written, 
in terms of the point process (3.3), as 

Sk,n = ^2 9n{X h V n )e R -i n x.({x- e 0^ d ) k : IMI > 1 , * = 1 ,- ■ ■ ,&}) ■ 

ieX \V n \,k 

Appealing to Theorem 3.2, 

Sk,n => Poi (V fc {x € (ML d ) k : ||x;|| >1, i = 1, ..., k}^j 



If we can now show that Lk >n -A 0 in probability, then the previous line suffices to 
prove the theorem. 

To this end, note that, in view of Theorem 3.2, the point process 
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converges to a non-trivial Poisson random measure. Also, it is easy to see from 
(3.4) that Rk, n /Rk+i,n ~> oo as n —> oo. Therefore, replacing the scaling con¬ 
stants R k+ i >n in (5.1) by R kn , the corresponding point process converges to zero 
in probability. Thus L k , n —> 0 follows, as required. □ 


e -kp-c Ei=1 <e.w> fc(o, y) dpj(6)d6dy 

S d_1 Jp>0, p+c^ 1 (6,yi)> 0, i=l,...,fc—1 y 


THEOREM 5.2. Assume the conditions of Theorem 4.1. 

(i) If a(Rk,n)/ r n —> c £ (0, oo] as n —> oo, then 

Pk-2{c(V n (1 B(0] Rk, n ) c ,r n )^J => 

p«(u I f 

7(]Rd)fc-l Jsd-1 Jp 

(ii) If a{Rk, n )/r n —> 0 as n —> oo, then 

Pk-2{c(V n nB(0-R k7n ) c ,r n )) 4 0. 

PROOF. For the proof of (i), note first that Theorem 4.1 (i). implies that 
(5.2) J2 g n (X i ,V n )e Xi(n) ({[0,™] d ) k ) ^ Poz (^ fc {([0,cx)] d ) fc }Y 

i£X \Vn\,k 

Due to the symmetry of the integral with respect to 6 € 5 d_1 , 
z^{([0,oo] d ) fc } 


1 

k\ 


- f [ [ e kp c 1 ^-‘ i = 1< ' 6 ’ yi ' > h(0,y) dpJ(0)d9dy 

' J(K d ) fc - 1 Jey 0 Jp> 0, P+C~ 1 {9,y i )> 0, i=l,...,fc—1 


2 d k\ 


/ 

Js d - 1 


e kp c 1 Sfer /r.(0, y) dpj{9)d6dy. 


p>0,p+c 1 (0,Ui)>0 , i=l,...,fc—1 


Further, due to the spherical symmetry of the density, the left-hand side in (5.2) 
has the same weak limit as 2 ~ d Sk, n - Consequently, S/c,n weakly converges to 


Poi ( j - f [ f e kp c 1 Si=i /j(o, y) dpJ(6)d9dy 

\k- J(m d )k-i Jsd- 1 Jp> 0, p+c~ 1 {0,y i )> 0, i=l,...,fc-l 

To complete the proof of (i) it now suffices to verify that L\. n -a 0 in probability 
as n —> oo, which follows along the same lines as in the proof of Theorem 5.1. 

To show (ii), note that 

4-2 [c(V n n B{ 0; R k ,n) c , r n )^j 
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< y ] h n (Xi) l( ||Xij || > Rk,n ) j — 1) • ■ ■ ) k ) + Lk,n . 

i£l \Vn\,k 

The above expression has the same weak limit as 

2 d M*i)<Uf)(([0,oo] d ) fc ) +L fc , n . 

i£X \Vn\,k 

Once again, L^ n —> 0 in probability, and so the assertion follows from Theorem 
4.i (a). □ 

5.2. Limit Theorems for Partial Maxima. We now turn to describing the limit¬ 
ing behaviour of the maximum distance from the origin of the random points con¬ 
stituting the geometric objects that we have been studying so far. More specifically, 
we consider the ‘maxima process’ 

(5.3) \/ l|JQlll ~ 4 ’ n , t € [0,1], 

Cl [* IY-) 

l<ll<-<i k <\Pn\t , 

9n (T .'Pn) = 1 

where a V b = max{a, b} for a, b € M, and Ck } n > 0 and c4 n £ I are the 
normalising sequences of the previous sections. 

For a concrete example, suppose that h n : (M. d ) k —y {0,1} is defined as in (2.4), 
where T is a connected graph with k vertices. Then, a A:-tuple of random points can 
contribute to the maxima process, only if its points serve as the vertices of graph 
isomorphic to T. 

Note that the process (5.3) only requires the computation of the maximum of 
|| Xi x || (suitably scaled and centred). Flowever, in view of (2.7), all the components 
in X\ = (X t] ,..., Xi k ) must be close to each other. Therefore, the results below 
will be robust as to which component is chosen from X- x . 

We start with the regularly varying tail case, which is essentially a corollary of 
Theorem 3.2. The limit in this case is a time-scaled extremal Frechet process (cf. 
[20]). The main difference between this and a classical extremal Frechet process 
is that the former exhibits dependence in the max-increments, while the latter does 
not. 

THEOREM 5.3. Under the assumptions of Theorem 3.2, let (ji,si) represent 
the points of a Poisson random measure with mean measure v\. x A/,., where if,- (.4) = 
Ok(A x (R d )^ _1 ) for a measurable set A C R d (is/c is defined in (3.5)), and 
Xk(B) = Afc({(zi,... ,Zk) G Lk , Zk € B}) for measurable B C [0,1]. Then, 

V <»-D[0,l], 

l<ii<~<i k <\Pn\t, k,n S[<t 

gn{XuPn)= 1 
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where D\ 0,1] is the space of right continuous functions from [0,1] into R with left 
limits. The limiting process is a time-scaled extremal (ak — d)-Frecliet process with 
finite-dimensional laws determined as follows: for 0 = Iq < t\ <■■■<! k < I, 
m > 0, i = 1,..., K, 

(5.4) 


V 

Sl<ti 

= exp 


< m ; i = l,...,K 


Sd-1 

{k\) 2 (ak - d) 7( R d)fc- 


K 


MO, y)dy ~ fk i- 1) ( A dj 

i<j<K 


i =1 


— (ak—d) 


PROOF. Restricting the domain of the point process convergence in Theorem 
3.2, we have 

X] 9n{Xi,V n )e^ x ^ (•) 

lGl \Vn\,k 

^X e 0o*d(-) inM p (([-oo,oo] d \{0}) x[0,l]). 

i 


The functional T : M p (([—oo, oo] d \{0}) x [0,1]) -A D[0, 1] defined by e^ Zun 
V n <- \\ z l ||> is almost surely continuous, (cf. p214 of [22].) Applying the continuous 
mapping theorem immediately gives the required weak convergence. 

What remains is to establish the precise form of the limit process, as in (5.4). To 
show this, note first that 

p( V IliHI <m,i = !,■■• ,k) 

Sl<ti 

( K ) 

= expi - G : ||z|| > /\ Vj}^) h{{U-i,U]) > ■ 

^ *=1 i<j<K ' 

By (3.5) and (3.6) we have that 


n k (\z <E R d : ||z|| > /\ rjj 

i<j<K 

$d— 


k\(ak — 


Ad f Mo,y)rfy ( A 

- d) J { Rd)*-1 V' ' 


Vj 


i<j<K 


— (ak—d) 


and 

from which (5.4) now follows. □ 
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We now turn to the case of light-tailed densities. Note firstly that a similar, but 
simpler, argument than the proof of Theorems 2.1 and 4.1 shows that under the 
conditions in Theorem 4.1 (i), the point process 


^2 9n{Xi,V n )e 

ieI \Vn\,k 


afRfc.n ) -1 (ll Xij ,j=l,...,k) , i/\Vn\ 


converges weakly to a Poisson random measure with mean measure /if- x A/, : . As 
usual, Xf. is the /c-dimensional Lebesgue measure concentrated on L}., and is 
given by 


Abe(Oo,bo] x ■■■ x (a fc _i,fr fc _i]) 

= tt [ I I l(aj < p + c _1 (6>,?/j) < bi, i = 1,..., k - l) 

k\ J {R d )k ~i j Sdl j ao 


k -1 

x expj -kp - c _1 ^2(9, yi) } h( 0,y) dpJ(9)d6dy 1 
1=1 


where a*, bi, i = 0,..., k — 1 are one-dimensional real vectors with —oo < a, < 
bi < oo, and J{9) = |ckc/c)#| is the Jacobian. Exploiting this result and mimicking 
the proof of Theorem 5.3, one can prove the following result, in which the limit is 
well described as a time-scaled extremal Gumbel process. 


THEOREM 5.4. Assume the conditions in Theorem 4.1, and let ( ji , sj) be the 
points of a Poisson random measure with mean measure p r X A/- , where fik(A) = 
Pk(A x (R^)^ 1 ) for a measurable set A C R and Xk(B) = Afc({(zi,..., Zk ) £ 
Lk , Zk £ J3}) for measurable B C [0,1]. 

{i) If a(Rk,n)/ r n ~> c £ (0, oo] as n —> oo, then 


(5.5) 


V 


Xii || Rk,\ 


l<ii<---<ik<\Vn\t , 
g n (Xi,Vn) = l 




V ii inD [o, l], 

Sl<t 


The finite dimensional laws of the limiting process are as follows: for 0 = to < 

ti < ■ ■ ■ < tx < 1. th £ R, i = 1, ■ ■ ■, K, 

JP > ( V 3 i<Vi, i = 

si<U 

= exp|- 1 / [ e~ c 1 ^=^ {9 ’ yi) h(0,y)J(9)d9dy 

{ («!) * J{ K d ) fc - 1 

K 

x - t£_i) e~ k ^<t< K % 

Z=1 
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(ii) If a(Rk,n)/r n —> 0 as n —>• oo, then the left-hand side in (5.5) converges to 0 
in probability. 

5.3. Limit Theorems for Partial Sums. Continuing on from the previous sub¬ 
section, we now consider limit theorems on partial sums. In particular, we shall 
focus on a stable limit case. As seen in a variety of related studies such as [12] 
and [13], proving stable limit theorems via a point process approach has become 
the gold standard. In order to obtain stable limits, however, the underlying random 
variables constructing partial sums must have infinite second moments. For this 
reason, we shall assume a regularly varying tail for the underlying density, and 
further, that the homogeneity exponent ak — d in (3.6) lies in the interval (0, 2). 
Combining this constraint with a > d and k > 2, we need to treat only the case 
1 < a < 1.5, k = 2, and d = 1. 

THEOREM 5.5. Under the conditions of Theorem 3.2, assume that 1 < a < 
1.5, k = 2, and d = 1. Suppose additionally that an indicator h : R 2 —>■ {0,1} is 
not only shift invariant as in (2.6) but also symmetric in the sense that 

(5.6) h(xi,X 2 ) = h(—x i, —xf) for all X\,X 2 G M . 

Then, Rfh |, 9n{%h 'Pn)Xj 1 converges weakly to a symmetric (2a — 1)- 

stable law. 

PROOF. Restricting the domain of point process convergence shown in Theo¬ 
rem 3.2, we find that in the space M p ([—oo, oo] \ {0}), 



converges weakly to a Poisson random measure with intensity if, where if (A) = 
f2 (A x M) for measurable dci 

Under our parameter restrictions, the homogeneity exponent in (3.6) is 2a — 1, 
and thus the limiting Poisson random measure can be represented in law by 


OO 



where ( rj , j > 1) is a sequence of i.i.d. Rademacher random variables taking +1 
and —1 with probability 1 /2, Tj is the j-th jump time of a unit rate Poisson process, 
and 
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Here, (rj) and (Fj) are taken to be independent. Notice that due to its symmetry, 

C a | rjTj <2r ' 1 j converges almost surely and has a symmetric (2 a — In¬ 
stable distribution. For more information about series representation of stable laws, 
see Section 1.4 of [25]. It is known that, for every <5 > 0, the functional Tg : 
M p ([—oo, oo] \ {0}) —>■ M defined by 

=J 2 zi1 (\ zi \ > 5 ) 

i i 

is almost surely continuous, (cf. Section 7.2.3 in [23].) Applying the continuous 
mapping theorem, we have, as n — > oo, 

R2,n J2 9n(Xi,'Pn)X il l{\X il \>R 2 , n 6) 

' 1&X \V n \,2 

OO 

=*► c a '£ r j rj < 2 *-»- 1 i ( C ' Q r i " (2 “ _1) “ 1 > 5 ). 

3 = 1 


As 5 4 , 0, we have 

OO 

(2a—I)” 1 


3 = 1 


t—( 2a—I) -1 


i(c Q r: 

Hence, it remains to show that, for every 77 > 0, 


> 5 )^c a Y J r j r; 

3 = 1 


—(2a—1)“ 


lim lim sup 1 

<540 n—yea 


J2 9n (Xi , V n )x h 1 ( 1X^1 < R 2 ,nS ) 

ieX |Pn|,2 


> r]R 2 



= 0 


However, by the Cauchy-Schwarz inequality, this will follow immediately if we 
can show that 

(5.7) lim lim sup e| V' g n {Xu 'Pn)X il l( \X h \ < R 2jn 5) ) =0. 

<540 71—>-00 ’ t. r— J 

ieI |Pnl,2 


We can write 

R 2 2 n^{ J2 5n(A' i ,F>„)X il l(|X n | <i? 2 , n 5)} 2 

{eI \Vn\,2 

oo 2 

= R 2 2 n ^P{|K|=m}E{ ]T g n {Xi,X m )X il l(\X il \<R 2>n 6)} , 

in =2 i £ 2 to ,2 


where = {A 1 ,... , A m } for m > 2. 
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Now introduce a triangular array of i.i.d. Rademacher variables (n, 1 < i\ < 
i 2 < oo), which are independent of (X,, i > 1). Then, by virtue of the symmetry 
of the Xi, it follows that, for all n > 1 and m > 2, 

(dn(Xii Xm) l(|^n | < , i ^ ^m,kj 

= (snOl^i, r i T’ m )r i Xj 1 l(|X il | < R 2 , n 5), i €lrn,k^j ■ 


We now observe that 

E{ £ 5n (A’ i ,A’ m )X il l(|X il | <f? 2)n ^)} 2 
i£X m) 2 

= ^ SVi(n^i , n^m) n^ii 1 ( \X U \ < R 2 , n S) | 

iEZm,2 

- £ £ E { Xi 1 Xj 1 l( \X^\ < R 2>n 6 , I-Xjj < 

iGlm,2jGXm,2 ^ 


X E^r;rj ■ L;fh m ) , rj<T m ) Tmj’ 


If i / j, n and ? j are independent, and so 


E 


■j^iTj j ^iX m ) fj n (‘^j ■ 1 jX m ) X m 


= E 
= 4 
= 0 


|^i 9n(^\Xi j Li2T m ) Tmj- E-^fj p n (rjTj , /'j<T m j <T m j- 

_1 , ^m) - 9n{-Xh ~XrrS) (sri(*j , X m ) ~ 9n(~X j, “^m)) 


where the last equality follows from (5.6). 

This implies that all cross terms will vanish, and so 

Ej J2 9n(Xi,X m )X il l(\X il \ <R 2 ,nS)} 2 
iGXm,2 

= Y ^{an{nx i , nX m ) x? i( \x h | < R 2 , n 5 )} 

iGXm,2 

= Y ^{an(X i, Xm) Xl 1( \X tl I < R 2 , n S) } 

1 dz^Lm , 2 

< QE{h n (Xi,X 2 )X 1 2 l(|X 1 | < f?2,n<5)} . 
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Finally, we have that 


R- 2 fn, XT 9n(X\,'P n ) < i?2,rU 

iSXi-p„ 1.2 


\V n \,2 
oo 


m=2 


< R^ n £ F{\V n \ = m) ("M E [h n (X 1 ,X 2 )Xf 1 (| X\ \ < R 2 , n 6)} 


n 


™\n - 

^2^, p R2,n& 


h n {xi,x 2 )xj l(|xi| < R2,n$)f{xi)f(x 2 )dxidx 2 


p2 

Il 2,n 




~ i? 2 


2,71 


/* -Ft2, n O r 

/ x 2 f(x) f(x\l + r n y/x\)h(0,y)dyda 
Jo J ]R 

r-R2, n S r f(x\l + r m y/x\) , 

/ x 2 f(x) 2 / sup-j^r-; - -h(0,y)dydx. 

JO J R m>l 


Here, the equality in the fourth line is obtained by the change of variables aq -o- x, 
x 2 -o- x + r n y. Since (r n ) is a bounded sequence, it follows from the uniform 
convergence of regularly varying functions of negative exponent (cf. Proposition 
2.4 in [23]) that 


sup 

m> 1 


f{x |1 + r m y/x |) 
f(x) 


1 > 


as s A 00 , 


uniformly in y G R with h{ 0, y) = 1. 
Thus, as x -A cx), 


U(x)= f sup + r™ y / x \) h( 0 ,y)dy ->■ [ h( 0 , y)dy G (0, 00 ) 
Jlm>i /(£) Jr 


This fact implies that x 2 f(x) 2 U(x) is regularly varying with exponent 2 — 2a. By 
an application of Karamata’s theorem (e.g. Theorem 2.1 in [23]), 


n 2 r n 
R'2.n JO 


rR2,nS 

/ x 2 f{x) 2 U{x)dx ~ ^-^-R2,nS d f{R2,n6) 2 U(R 2 , n S) 
Jo 3-2 a 


n 2 r n 


n 2 r n R2,nf(R2 
d 3 ~ 2a 


S 3 


-2a 


[ h(0, y)dy 
Jr 


3 — 2 a 


[ h(0,y)dy 
Jr 


3 — 2 a 

as n —j cx). 


Since <5 3 2 “ —>• 0 as <5 4- 0, (5.7) follows. 


□ 


6. Appendix, 
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6.1. Proof of Theorem 2.1. For the completion of Theorem 2.1, we need sev¬ 
eral important ingredients, all of which belong to the ‘Palm theory’ of Poisson 
point processes. They are useful when computing expectations related to Poisson 
point processes. We recall abbreviations (2.1), (2.2), together with (2.3). 

LEMMA 6.1. (Palm theory for Poisson point processes, [2]; see also Theorem 
1.6 in [21]) Let (X,, i > 1) he i.i.d. W 1 -valued random variables with common 
density f, and for n > 1, let V n denote a Poisson point process with intensity nf. 
Let u n (y ,X) be a measurable function defined for y = (yi,, yf), y% £ and 
a finite subset X D f of d-dimensional real vectors. Then, 



where V' n is an independent copy ofV n . 

LEMMA 6.2. Under the same notation as Lemma 6.1, let u n (y) be a measur¬ 
able function defined for y = (yi ,..., yk), yi £ R d . Let be the k-dimensional 
Lebesgue measure concentrated on the upper diagonal part of the unit cube; viz. 
Lk = { ( 21 , ..., Zk ) £ [0, l] fc : 0 < z\ < ■ ■ ■ < Zk < l}. Then, for a measurable 
set B C Lk, 


E { XT u n{Xi)l(}/\Vn\ £ B)\ ~ ^E{u n (X k )} X k (B), flin AOO. 


Proof. 




00 


m=k 


The proof then follows from the fact that 



□ 
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LEMMA 6.3. Under the same notation as Lemma 6.1, let u n (y) be a non¬ 
negative measurable function defined for all finite subsets y = (y \...., y k ), yi G 
W l . Suppose further that n k ¥,u n (X k ) —> C, n —> oo for some constant C > 0. 
Then, for a measurable set B C Lk, 


E 


Y u *(*0 G B) - Y MM l(i/n G B) 

i€l\-Pnl,k i£Zn,k 

Proof. We can proceed as follows. 


0 , n —> oo . 


E Y MM i^/l^nl € B) - Y M* i) M/n G B) 

oo 

= Yj P{l^ re l = m } #{i G n mB} - #{i G n nB} E {u n {X k )} 

m= 0 

oo 

~ Cn~ k Y ^{\^n\ = m} #{i € T n)k n nB} - #{i G T mM n mB] 

m =0 

OO fc 

< Cn~ k Y HIM = rn} | #{i G Z n>fc n nB} - k (B) 

m= 0 

oo x ( 

+ Cn- fc ^P{|P„| = m} 


n k — m k 


MB) 


m =0 
oo 


fc! 


+ C7r fc P{|P ra | = m} — Afe(S) - #{i G x m)fc nmB} 

m =0 

= Ji + Ji + J3. 


Evidently, J\ —> 0 as n —> 00 . It is also easy to check that 


j 2 = c^e 

k\ 


\Pn 


n 


- 1 


0 as n —> 00 . 


As for J3, given e > 0, there is an integer N > 1 such that 


fc! 


—t #{i G n mB} - \ k (B) 


m 


< e for all m > N . 


Then, there exists a constant C' > C such that for all m > N, 


N -1 


J 3 < C'n~ k Y P {\M = m } r -JT + Cen ~ k Y = 


m =0 


m=N 


fc! 
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Obviously, the first term on the right-hand side vanishes. The second term is bounded 
by Ce(k\)~ 1 n~ k E\V n \ k . Since sup n>1 n~ k E\V n \ k < oo and e is arbitrary, we con¬ 
clude that J\ + J 2 + J 3 —> 0 as n —> 00 . □ 


Before starting the proof of Theorem 2.1, we note that we shall often use, with¬ 
out further comment, the fact that, for every k > 1 , 

f n\ n k 

\k) ~ IT’ 

in the sense that the ratio of the two sides tends to 1 as n —> 00 . 

Proof of Theorem 2.1. For the convergence of N^\ according to Kallen- 
berg’s theorem (see Proposition 3.22 in [22]), it suffices to show that 


(6.1) 

E {NW(R)} 

-> E{iV (fc) (i2)}, 

and 



(6.2) 

P [n^(R) = 0} 

-A P jiV (fc) (.R) = 0} 


for every disjoint union of measurable sets of the form R = Up=i( j 4 p x B p ), 
where A p is relatively compact in E). with 1/^(8A p ) = 0 (i.e. the boundary of A p 
has ^-measure 0), and B p is a measurable set in L&. 

For the proof of (6.1), we can set, without loss of generality, m = 1, and write 
A = Ai, B = B 1 . By virtue of Lemma 6.2 in the Appendix, together with (2.17), 

E{Ni k \A x B)} ~ {h n (X k ) = 1, £ A }\ k (B) 

^ MA)Xk(B) =E{N^ k \A x B)} . 

Now we proceed to (6.2). Setting 


£i ,n — h n (X i) 6 ^^.(n) \ (-^) ' * — {ft, . . . , G ^-n,k : 


we see that 


F{N^(R) = 0} - ¥{N^ k \R) = 0} 

< F{N^(R) =0}-F{ ]T ^i,n = 0} 

+ |P{ E £i,„ = 0} - ¥{N^(R) = 0} 
i £Z n ,k 


= I\+h. 
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We recall the following basic inequality: for integer-valued random variables X 
and Y defined on the same probability space, 

|P{X = 0} - P{Y = 0}| < E\X - Y | . 

Combining this inequality and Lemma 6.3 in the Appendix proves I\ —> 0 as 
n —> oo. To demonstrate that I 2 —> 0 as n -A 00 , we introduce the total variation 
distance: for real-valued random variables X and Y defined on the same probabil¬ 
ity space (fl, J 7 , P), 

d TW {X, Y) = sup \P{X € A} - P {Y G A}\. 


Using this norm, 


(6.3) 


P{ ^ Ci,n = 0}-nN {k \R) = 0} 


< dxv 





Since the total variation distance between two Poisson variables can be bounded 
by the difference of their means (see Lemma 7.3 in [8]), we have 

<kv(Poi( e{ J2 6,n}), < e{ . 

LT n ,k 

Thus the final term in (6.3) converges to zero as —>• 00 . 

To handle the first term on the right hand side of (6.3), note first that each £i >n is 
a Bernoulli random variable. To handle dependencies, create a graph with vertices 
the indices i G 'I n .k by placing an edge between i and j (write i ~ j) if i and j share 
at least one component (i.e. if |i n j| > 0). Then, ~) provides a dependency 
graph with respect to (£ii G that is, for any two disjoint subsets I\, 

I 2 of T n ,k with no edges connecting I\ and 1 2 , (hi.n, i G / 1 ) is independent of 
(^i n , i G 12 )- Therefore, we are able to apply the so-called Poisson approximation 
theorem (see [2] and also, Theorem 2.1 in [21]): 

dxv (E £i,n, 6,n}) 

iel n ,k 
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< 3 l E E E Ui,n}E{4n} + E E E te,-4n} \ > 

^ iG^n.fc jgA; iGX n ^. jsNi\{i} 

where Ni = {j € I n ,k '■ |i n j| > 0}. Now, for large enough n, 


E Ui,n} < E p {^(^) = 1, 4 n) e A p } < 2 

P= 1 


— 1 m 


El ^k^Ap), 


P= 1 


and thus 



—2 m 


(E Uk ^l 


P= 1 


p=l p'=1 


E { 6 ,nCj,n} < E E P \ M*fc) = E ^ € ^ 4 pj 


^n(Xi, • • • j A'/jXfc+i, . . . , X 2 k-l) — 1, 

( c fc,n( x i “ d k}n S(Xj )), i = 1, • • •, l, k + 1,..., 2k 



Regardless of the definition of E k , there exists a compact set K\ C E-ik-i such 
that for all p, p' = 1,..., m. 


{(xi ,...,x 2 k-i) € (R d ) 2k 1 : (x 1 ,...,x k ) G A p , 

(x\, ,Xi, Xk-\-\i • • • 5 %2 k—l ) ^ -^1 


Thus, for some constant C\ > 0, 

El El E {£i,ra£j,n} 

jeN-Mi} 

^{^.nlidinjl = 0 

fc-i 

< CiEn 2H P{hn(^) = 1, M*t, • • ... ,X 2fe _0 = 1, 

(=i 

X 2 k-i £ Xi | —> 0 as n —> oo , 
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where the last convergence follows from (2.18) and now, (6.2) is proved as re¬ 
quired. 

~(k) (k) 

In order to prove that Nn has the same weak limit as Nn , we only have to 
verify that for every non-negative continuous function / : E k x L k —v K + with 
compact support, 

Ni k) U) - «?>(/) 

= ( h n(^)~Sn(Xi,V n ))f{Xl n> , \/\P n \) 40 . 

i&x \r n \,k 


Let K 2 be a compact set in E k so that the support of / is contained in K 2 x L k . 
Noting that [|/||oo = sup (a . |W)eXax£fc f(x,y) < oo, 

< ll/lloo l(hn(Xi) = 1, g n (Xi,P n ) = 0, X^ EK 2 y 

i£l \v n \,k 

In view of the Palm theory in Lemma 6.1 in the Appendix, we need to show that 

n k F{h n (X k ) = 1, 9 n{X k , X k U V n ) = 0, Xl n) G K 2 } ^ 0 , 

where V' n is an independent copy of V n . Combining (2.17) and (2.19) completes 
the proof. □ 

6.2. Proof of Theorem 3.2. 

Proof of Theorem 3.2. By Theorem 2.1, we only establish that the three 
convergence conditions, (2.17), (2.18), and (2.19), with c k n , d k)Tl replaced by Rk. n 
and 0, respectively, are satisfied. Since the proof of (2.17) is very similar to (and 
actually even easier than) that of (2.19), we only check (2.18) and (2.19). We shall 
start with (2.19). In view of the Portmanteau theorem for vague convergence (e.g. 
Proposition 3.12 in [22]) we need to show 

(6.4) Q P{ g n (X k , X k U V' n ) = 1, Rf]X k € K } 4 u k {K) 

for all relatively compact K in E k = ([— oo, oo] rf ) fc \ {0} for which v k {dK) = 0. 
Without loss of generality, we can take K = (oo, f>o] x • • • x {a k - 1 , b k - 1 ] C E k , 
where o,, bi, i = 0,..., k — 1 are d-dimensional real vectors. Recall that any 
relatively compact set in E k is bounded away from the origin, and as assumed in 
(2.7), h(x i,..., x k ) = 1 only when aq,..., x k are all close enough to each other. 
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Therefore, we can and shall assume that each (a,;, bj] is bounded away from the 
origin. In particular, we assume that there exists an p > 0 such that 

(6.5) (dj, bj\ C {x G W l : ||x|| > 77 } , i = 0, ..., k — 1. 

Consequently, we have 


P{ 9 n{Xk, X k U V' n ) = 1 . dj_i X R k ] n Xi A 6 j_i, i = 1 ,..., k } 
= E^h n (X k ) l( Oi_i -< A &»_i, i = 1 ,... 

xP G(X)., r n ) is an isolated component of G(X k U V ' n , r n j 



hn{xu ■ ■ ■ ,Xk) l(oi_i A A bi- 1 , 7 = 


x exp{—np(xi, ...,Xk] r n )}f(x 1 )... f(x k )dx. 

Let //. denote the last integral. The change of variables x\ ftx + r„y,_ 1 , 

i = 2,... ,k, together with the location invariance of h, yields 


4 = 


(l) r n k 1] [ . /(*) / ^(0,y) n l(oi A 4 fc ^(x + r n ?/i) A 6i) 

W Jao^R-^bo J(Rd)k-l 

X /(x + r n yj) exp{—np(x,x + r n y; r n )}dydx. 


Applying the polar coordinate transform x -H- (r, 0) with .7(0) = (7.x/47|, and 
changing variables by setting p = r/R ktfl , gives 


4 = Q r n' k ~ 1 ^k,nf (4fc,n e t) fc 

/* /* /*oo 

x / / / l(a 0 A A 6 0 )p^ 1 /(4fc, n ei) _1 /(4fc,npei) 

•7(R d ) , =- 1 JS' 1 - 1 do 
fc -1 

x 1 (a* A + r n yi/R k!n A &*) /(i? A . ! „ei) _1 /(4 fc>n ||p6» + r n yi/4 fcin ||ei) 
2=1 

x exp{-np(R ktn p9, Rk, n p0 + r n y; r n )} h(0,y)dpJ(9)d9dy, 

where J{9) = \dx/d6\ is the usual Jacobian and S d ~ l denotes the (d — 1)- 
dimensional unit sphere in M d . Recalling (3.4) and the assumption that / has a 
regularly varying tail of exponent —a, we see that for all p > 0, 9 £ S d 1 and 
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( 6 . 6 ) 



k -1 


P@ “ 1 “ ^nVi/Rk,n — 


X 


f (-^/c,n||P^ “1“ ^ nVij -^/c,n||^l) 


^ fc—1 

- . 



2 — 1 


Subsequently, we shall show that 


(6.7) 


np{Rk, n pO , -Rfc.nP# + f n y; r n ) -> 0, n ->oo 


for every p > 0, 8 € S d 1 and y,; e M d , i = 1,..., fc — 1. By the change of 
variables, 

( 6 . 8 ) 

np(R k , n p8 , Rk,nP@ + f n y; T n ) 



Appealing to the Potter bound (e.g. Theorem 1.5.6 in [7]), for every £ € (0. a), 
there is C\ > 0 such that for sufficiently large n, 


f{Rk,n\\pO + r n z/Rk, n \\ei) 



zGB(0; 2)UUi4 1 1 B(yi\ 2) ' 

< Ci(2p-(“ + ^ + 2p-( Q -^) . 

Therefore, the supremum over n of the integral in (6.8) is finite. On the other hand, 
(3.4) ensures that nr d f(R ktn e 1 ) —> 0 as n —> oo and hence, the convergence in 
(6.7) follows. 


Note that 



fc-i 


P 


A—l—ak 



2—0 
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If 

W- J( R d ) 


h(0,y)dy / ||x 

Jdi^x^bi i=0,...,k —1 

= v k ((a 0 ,b 0 ] x ''' x (a fc _i,6fc-i]) • 


I — ak 


dx 


Thus, the proof of (6.4) can be finished, provided that the convergence in (6.6) 
holds under the integral sign. Indeed, one more application of the Potter bound, 
together with (6.5), verifies the following: for every £ S (0, a — d), there exist C 2 , 
C 3 > 0 such that 


l(a 0 < p9 < b 0 ) 


f(Rk, n pe 1 ) 
f(Rk, n ei) 


<C 2 l(p>p) (p-^+O+p-M), 


and 


k -1 

^ pd + r n yi/Rk : n C hi) 

2—1 


f (Rk,n\\p9 C Tnl/i/Rk,n\\&l) 

f(Rk, n e 1 ) 




Since f^°(p d 1 + p d 1 “ £)dp < 00 , the dominated convergence theorem 

justifies the convergence under the integral sign, and so we have completed the first 
part of the proof; viz. (2.19) is satisfied. 

Next, we shall prove (2.18). For l = 1,..., k — 1 and every compact set K C 
E 2k -i = ([— 00 , 00 ]“) \ {0}, we can assume without loss of generality that 

there exists rf > 0 such that 

K C {(xi,.. .,x 2k -i) G (R d ) 2k ~ l : ||z,|| >rf, i = 1,... ,2k - 1} . 

Proceeding in the same manner as above, the probability in (2.18) is bounded by 

n 2k ~ l v{h n (X k ) = 1, hniXu .. .,X h X k+1 ,.. .,X 2k _ t ) = 1, 

||2Q|| > rj'Rk , n , i = 1,. • ■ ,2fc - l j 

= n 2 k - l rd(2 k - l -i) R d n f { R knei) 2k- 1 
P P P OO 

X / / / p d ~ l f{Rk, n ei)~ l f{Rk,nPe\) 

J(j & d)2k-l-l Jgd -1 
2 k-l-l 

1 (||p$ 4" ^nyi/Rk,n [| ^ V f (Rk ,n& l) /(^fe,n[|P^ 4“ ^n2/i/-^fe,n||si) 

i=l 

x /i(0, ?/i,, 2/fc-i)/i(0, 2 / 1 ,..., 2/z-i) 2/fcj • • ■, V 2 k-i-i) dpJ(9)d9dy . 


The asymptotic order of the last integral is O (n 2k l r d ^ 2k 1 ^ R d n f{R kn e i) 2fc ; ), 
which tends to 0 as n —>• oo for every / = 1..... A: — 1, and so we are done. □ 
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6.3. Proof of Theorem 4.1. For the proof of Theorem 4.1, we need a prelimi¬ 
nary lemma. 

LEMMA 6.4. (Lemma 5.2 in [3]) Given (Rk,n, n > 1) as in (4.5), let ( q m (n ), m > 

0, n > 1) he defined by 

4m(jt) — 0‘{,Rk,n) iffi (pP{Rk,n) T tTlj Rk,n^j ; 

equivalently, 

fi{Rk,n + a(R k ,n)qm(n )) = f>(R k ,n ) + m 

(since ip is an increasing function, the inverse function ip*~ is well-defined every¬ 
where). Then, given e > 0, there is an integer N e > 1 such that 

<lm(p) < e m£ /e for all n > N e , m> 0 . 

Proof of Theorem 4.1 . The proof is rather long, and so we break into two 
main units, one each for the two main cases, further dividing the proof of the first 
case into three sign-posted parts. Hopefully this will help the reader to navigate the 
next few pages. 

Proof of statement (z) 

Part 1 : As argued in the proof of Theorem 3.2, we need only check (2.18) and 
(2.19). As for (2.19), we need to show that 

(6.9) (fj P{ g n {X kl X k U V n ) = 1, X^ G K } -A u k (K) 

for all relatively compact K in E k = ((— oo,oo] d ) A for which v k {dK) = 0. 
Without loss of generality, it suffices to check the case K = /q x • • • x I k _ i, where 
Ii = (a t , //,] c Ml 1 and a, , b, are d-dimensional vectors such that —oo -< a, A A 
oo. Define 


H = S^ 1 n {(* 1 , ..., z d ) € R d : zi + • • • + z d > 0}, 

and let H c = S d " 1 \ H. Let ©i = ATi/||Ai||. Since a(R k , n )~ l R Kn S(X i) -»■ oo 
a.s. (in a componentwise sense), we have that, for large enough n, ©i G H c and 
a(R kt n)- 1 {X 1 - R kjU S(Xi)) y a 0 (y —oo) do not occur simultaneously. Thus, 
the left-hand side in (6.9) equals 



4uK) 


k 

1, ©! € H , xff 1 ^ G J] /j_! | + o(l) 

i=l 
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as n —> oo. Letting J k denote the leading term in the above, we can write 

k 


J k = (JJ 1( x k e n Ii- 1, ©i G H) 

x p|g(<L/ c , r n ) is an isolated component of G[X k U V' n , r n ) X k 

?) / h n (x\,..., x k ) 

K J J{R d ) k 

x l^a(f?fc,n) _1 (^ - Rk,nS(xi )) € /j -1 ,i = l,...,k, xi/||xi|| € 
x exp{—np(zi,... ,x k \r n )} f(x i)... f(x k )dx., 

where the definition of p was given in (3.7). 

The change of variables x\ 44 x and x t 44 x + r n y,_ |, i = 2, ... ,k yields 




_ n \ ^rf(fe-i) 


l(a(-Rfc,n) 1 ( a: - Rk, n S(x)) e J o , z/IM| € f(x) 

~ k —1 

x / M°,y) TT l( a(74 n ) _1 (x + - R k , n S(x + r n yi)) G /*) 

V^- 1 f=i v 

x f(x + r n yi ) exp {-np(x, x + r n y; r n )}dydx. 

Further calculation by the polar coordinate transform x 44 (r, 0) with J{0) = 
\dx/dO\ and the change of variable p = a(R ktn )~ l [r — R k , n ) gives 

Jk = (l) r% k -Va(R k , n )R d k ~ 1 f(.R k , n e 1 ) k 


k-i j S d-1 j 0 ^ 


"J Li h( 0, y) dpj(9)d9dy , 


where 


Lx = l( P e-a(R k ^)- 1 R k , n (S(e)-e) e I ih 9 gH) 


L 2 — ( 1 + 


d -1 


0/(Rk,n 

Rk,n 

-^3 f (RkjnP l) / ^ “ 1 “ Q , (fik,n)p) ^ 1 ^ 5 

fc-1 

La = n 1 ( P° + ~TW~v Jl 
11 V a(R ktn ) 


i=l 
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R kj 


^{.Rk,r, 


s{{Rk,n + a{Rk,n)p)0 + r n y^j -9 G I^j , 


fc-i 


l 5 = Y\f(Rk, n ei) 1 f [\\{Rk,n + a{Rk,n)p)0 + r n yi\\e^j , 

1=1 

Lq = exp j -np((R ktn + a(R ktn )p)9 , (R k . n + a(R k)n )p)9 + r n y; r n ) | 


Part 2: For the application of the dominated convergence theorem, one needs to 
compute the limit for each Li, i = 1 ,..., 6 , while also establishing finite upper 
bounds for each term. We shall begin with the indicators L\ and Ij\, which are 
trivially bounded. The fact that a(R ktn )~ 1 R k ,n oo ensures that 

L\ —y 1 ( p6 G Iq , $ ^ 0), rr —^ oo . 


Observe also that 


Li < l(p > —M') for some M' > 0 . 

To see this in more detail, choose M' > 0 such that mini<j<,/ a,^ > -M' (the 
superscript (j) denotes ‘jth component 1 of a given vector). Fix p and 9 such that 
Li = 1. Then for all j with 9^ > 0, 

P > p eU) = pQ(3) - - OW) > a { 0 j) > -M '. 


Note that 9 G H guarantees that at least one component in 9 must be positive. 

In what follows, we prove the assertion when M' = 0. The proof for a general 
M' is notationally more complicated, but essentially the same. 

Before moving to L 4 , note the following useful expansion, which will be applied 
repeatedly in what follows: For each i = 1,... ,d, 

( 6 . 10 ) 


( Rk,n T O j {.Rk,n)P) @ T TnUi — Rk,n T d(,Rk,n) f P T 


(Q,yi) + in{p,Q,yi) 

a {Rk,n)/ r n 


so that 7 n{p,9,yi) —> 0 uniformly in p > 0, 9 G S d ~ l , and \\y l \\ < M ( M is 
determined in (2.7)). 

Turning now to L 4 , we shall rewrite the expression within the indicator as fol¬ 
lows: 


( 6 . 11 ) 


p9 + 


(^{Rk,n) 


Vi + 


Oi r< 

Pn 


where 
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OLn = {{OiVi) +ln{p,0 ) y i ))6 + 


Rk.i 


-“7c,n 


— m 


with 

and 


m 


(*) = (I 




Rk,n 1 

| ^(-^/c,n) 

r n ' 

V -^/c,n 

(* W , 

..,x (d) ) ef 


p)0 +Vi 


Pn — 1 + 


-Rfc.n 


P + 


(0,Vi) +ln{p,0,yi) 


As seen above, L\ —> 1 (p9 G Iq, 9 A 0) as n —> oo, so it is enough to discuss 
the convergence of (6.11) for every p > 0, 9 A 0, and \\y L || < M. Then, for large 
enough n. 


Rk'i 


1 + 


. n 


p) o + yi y 0. 


Thus, we have that a n —> (9 , y,,)() — yi and j3 n —>• 1 as n -A oo. Now we have 

and in conclusion, for every 9 A 0 , 


fc-i 

L 4 —> J^[ 1 (^p + c 1 (6,yi))9 G ■ 

i=1 

Regarding L 2 , it is clear that for every p > 0, L 2 -A 1 as n -4 oo and it is also 
easy to check that L 2 < 2(p V l) d_1 . 

As for L 3 , we write 

L 3 — n T a(f?fc in )p) exp | 1p(Rk,n T ®(^fc,n)p) T • 

An elementary calculation (e.g. pl42 in [14]) shows that 

/A 1 a {Rk,n) 

(6.12) —- - -— -A 1 

^\Rk,n T &{Rk,n)'V) 

uniformly on bounded u-sets. Namely 1/a is flat for a. Therefore, for every p > 0, 

exp {-ip(Rk, n + a(R k ,n)p) + V’CRfc.n)} -A e~ p , n -A 00 . 
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Since L is flat for a, it follows that L 3 — > e~ p as n — > 00 for every p > 0. 
For the upper bound of L.> on {p > 0}, we apply Lemma 6.4. Choosing e £ 
( 0 , (d + 7 fe) _1 ) and recalling that ip is non-decreasing, 

exp| -'tp(R-k t n + a(Ri^n)p) + ip{R k ,n)} Hp > 0 ) 

OO 

= X: 1 (?rn(n) < P < g m+ l(n)) exp{-^(.Rfc,n + a(Rk,n)p) + 1p(Rk,n)} 
m=0 

OO 

< 1(0 < p < £- 1 e (m+1)e ) e - m 

m=0 

for all n> N e . 

On the other hand, using the bound in (4.3), on {p > 0}, we have, for sufficiently 
large n. 


L^Rk^) 1 L(R k>n + a(Rk, n )p) < C (l + ^nlpY < 2 C(p V 1 )U 

V / 

Multiplying these bounds together, we have 

OO 

L 3 l(/9 > 0) < 2C(p V iy !(0 < P < e- 1 e( m+1 ) £ ) e~ m . 

m =0 


Next, we turn to L 5 . First, denote 


(6.13) 

Since c = limn^oo 


€n(p,0,y) = 


{9, yi) +7 n{p,0,y) 


O j (Rk,n') /'iV; 
a(Rk } n)/ r n is strictly positive, 


(6.14) A= sup |^ n (p, 0,y)| < 00 . 

n>l, p> 0, 
eeS 11 - 1 , ||j/||<M 


Using the expansion ( 6 .10), we can write 

fc-t 


L 5 = n L ( R k,n) 1 L^R k>n + a(R k)n ) (p + f n (p, 0, yP)) ) 

rpHMyp a (i2 fc ) 


Z=1 


x exp < — 




dr > . 
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Due to the uniform convergence (6.12) and (6.14), for every p > 0, 9 <E S d 1 , and 

hi\\ < m. 


L 


P+£n(p,6,yi) 


Q’i.RkyT, 


d^Rkyn 4“ o(-Rfc,n)^) 


dr -A p + c l (9,yi), n —> oo 


and 

HRkyn^L^Rkyn + a(R k ,n)(p + Cn(p,9,yi))^j -a 1, n ^ oo. 

We thus conclude that 


fc-i 

L 5 -a exp{ —(/c - l)p - c' 1 ^(6 , ,y i )} 

2=1 


for every p > 0, 9 € S'** -1 , and ||y;|| < M, i = 1,..., k — 1. 

To provide an appropriate upper bound for L-, on {p > 0}, note that, for large 
enough n, 


k-1 

\ exp 


2=1 



P+€n(p,0,yi) 






dr > l(p > 0) 


< ex | V /° 

\ J—A &(-^fe,n 

fc_1 rP+^(Pi 0 >yi) 


-§i 

C e 2A(fe-l) 




“ 1 “ tz(-Rfc j 7 l )r) 


dr l(-A < p + £ n (p,0,yi) < 0) 


dr l{p + i n {p,d,yi) > 


It follows from (4.2) and (4.3) that there exists Ci > 1 such that, on {p > 0}, 

L^Rk^n) R (^Rk,n “1“ ^{Rk,n) (p “1“ £n(P? ^5 2/i)) ^ 

= i(^fc,n) _1 ^(^fc,n + a(^fc,n)(/0 + ?n(p,6 l ,yi))) 1 (/> + £n(P, Vi) > --4) 

< 2 + L(Rk t n ) 1 L^Rk ! n + a(Rkyn)(p + £n(p ; 0) 2/i))) l(p + £n(p ; Vi) > 0) 

<2 + c( 1 + ^d( P + f„(P,«,!/ i ) )) 7 

< 2 + 2C|l V (p + £n(p ; 0) Pi))) 

< 67f- 1)_1 (pvir 







POINT PROCESSES UNDER GEOMETRIC CONSTRAINTS 


49 


for sufficiently large n. This in turn implies that 


k -1 

H L{R Kn )~ l L[R Kn + a(R kjn ){p + UP, 0, Vi))) < Ci (p V 

i=1 

on {p > 0 }, and further, 


l 5 < Ci e 2j4 ( fc_1 ) ( p v iy( k -v 


holds on {p > 0 }. 

Finally, we turn to Lq. 

- log L 6 = np((R k , n + a{R k , n )p)0, ( Rk,n + a{Rk, n )p)0 + r n y; r n ) 

= nr%f(R kjn e i) [ f(Rk, n e i ) _1 

■/b( 0 ; 2 )uU^ %; 2 ) 

X / ^ + a(R k ,n) (P + £n(p, 0, z))) ei^j rfz ■ 

As argued in the derivation of the bound for L 5 , the supremum of integrand over 
n can be bounded by a constant multiple of (p V l) 7 . On the other hand, (4.5) 
guarantees nr^f {R k ^ n e{) —>• 0 as n —> 00 and hence, Lq —> 1 as n -4 00 . 

From the argument thus far, it follows that for every p > 0, 0 £ S' 1 ' 1 , and 

Hz/ill < M, * = 1,... ,d, 
f[ Li -+ e -k^- 1 'L£iHe,y i ) 

i— 1 

xl (pOelo, 6 y 0, (p + U(0,yi))9 <E Ii, i = 1,... ,k - l) . 

To finish the argument, it remains to check the L 1 -integrability on {p > 0} of the 
upper bound for nf=i U As we have seen so far, 

6 

n^>°) 

i=i 

00 

< 4CCi e 2j4 ( fc_1 ) ( p v l) d - 1+7fe 1(0 < p < £- 1 e( m+ 1 ) e )e _m 

m =0 


for sufficiently large n. Indeed, 



d— l+7/c 


^2 1(0 < p < e -l e (™+lF) e -m (ip 
m =0 
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< 


e (d+7fc)e 

e d+7fc 


e -[l-(d+7fc)e]m 

m=0 


so that the right-hand side is finite because we took 0 < e < (d + 'yk)~ 1 . Applying 
the dominated convergence theorem as well as (4.5), it turns out that J k —»• u k (K), 
n —>• oo as required. 

Part 3: Next, we shall prove that (2.18) is satisfied. For l = 1,..., k — 1, and any 
compact set AT C E 2k -i = ((—oo, oo] rf ) 2fc \ there exists B > 0 such that 

AT C {(xi,.. .,x 2 ; fc-0 6 (M d ) 2fc_z : xy A -SI, i = 1,..., 2A; - l } , 


where 1 is a d-dimensional vector with all entries 1. Then, the probability in (2.18) 
can be bounded by 

n 2 fc -'P{M* fe ) = 1 • h n (X u .. -,X h X k+1 ,.. .,X 2k _0 = 1 , y -SI} 
= n 2k ~ l / h n (xi,...,x k )h n (x 1 ,...,xi,x k+1 ,...,x 2 k-i) 

J(jg,d)2k-l 

X 1 (xi- R k , n S(xi) y -a(R k ,n)B 1, i = 1,... ,2k - l, xi/||xi|| € H^j 
x f(xi)... f{x 2 k-i)dx. + o(l), n —)■ oo . 


The equality above follows from the fact that || € IP and A| — R kn S[X\) > 

—a(Rk, n )Bl do not occur simultaneously. 

By precisely the same change of variables and polar coordinate transform as in the 
proof of (2.19), the leading term of the last line above equals 


n 2fc-i r d ( 2fc-I- 1)a(i2fcjn)i2 d- 1/(i2fcjnei) 


2k—l 


[ / / 

J( R d)2fc —i —i J S d -1 J o 

d—1 


1 (p0- R ^ n ' (S(6) - e) A-SI, 6£H 




X 1 + 


^{Rk^n) 
Rk.’i 


Pj f (Rk,n£ l) 1^ 


2k—l—l 


n i S (, + 


i =1 






Rk,r< 


^ ^ (Rk,n @ @ 

X f (-Rfc,n^l) f | (-Rfc,n “1“ tt(Rk,n)p) @ H - ^n2/z| 


>- -SI 


x /i(0, yi, ..., 2/fe-i) ^(0,2/i , — , yi-uVk, • • •, 2 / 2 ife—z—l) dpJ(6)d6dy . 










POINT PROCESSES UNDER GEOMETRIC CONSTRAINTS 


51 


The asymptotic order of the above expression is 

0(n 2k ^ l r ^ 2k ~ l ~ 1 ^a(Rk t n)R < k~n f(Rk,n ei ) 2k l ) j 

which vanishes as n —> oo for every l = 1 ,..., k — 1 , and so the proof of (i) is 
complete. 


Proof of statement ( ii ) 

We need only show that, as n -A oo, 

Ni k \f)= M*i)/(*i (n) , i/im)Ao 

ie2: |'Pr l |,fc 

for every continuous non-negative function f : Ek x L k M + with compact 
support. Note that the support of / is contained in 

{(xi,...,x fc ) e ( M. d ) k : Xi P -Bil, i = l,...,k} x L k 

for some If > 0. Therefore, 

Ni k \f)< ll/ll* X] l{K(X i ) = l,X^y-B 1 l). 

' l€X \Vn\,k 

Observe that X l;j — R kri S(X t . ) A —a(Rf. n ) If 1 implies 

|| Rk,n ^ Q‘{R k ,n)B 2 

for some B> > If - Hence, 

Ni k \f) <11/1100 £ = 

ll-X’ij II - -Rfc.n > -a(Rk,n)B 2 , j = 1, • • •, k) . 

In view of Lemma 6.1 in the Appendix, what now needs to be verified is that 

n k P{h n (X k ) = 1, ||Xj|| -R n > -a(R ktn )B 2 , i = 1,..., k } -A 0 . 


Once again, applying the same kind of change of variables and polar coordinate 
transform, together with (6.10), 

n k F{h n (X k ) = 1, ||Xj|| - R n > -a{R k ,n)B 2 , i = l,...,k} 


n k r d n ^a(R Kn )R d k ^f(R Kn e i) 



a(Rk,n) V - 1 
Rk,n / 
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k -1 

X f{Rk,ne-i)~ l f[[Rk,n + a(R kt n)p)e\) n 1 {p + €n(p,0,yi) > -B 2 ) 

2=1 

X fiRk.nei)- 1 f(jRk,n + a(Rk,n){p + Cn(p,d,yi))^ei \ h( 0, y) dpJ(9)d9dy, 

where the definition of £ n (p, 9, y L ) is given by (6.13). 

As argued before, 

f 1 + a(Rk^p\ d 1 < 2 {pyl) d-i 

\ -K'k,n J 

and, on {p > —B 2 }, 

L{Rk,n ) 1 L(Rk,n + 0>{Rk,n)p) 5; 2(7(p V l) 7 . 

Furthermore, since a is eventually non-increasing, 

expj -'ip(R k ,n + a(Rk,n)p) + 'ip(Rk, n )} <e~ p . 

As argued before, it follows from (4.2) and (4.3) that there exists C 2 > 0 such that 
on {p > -B 2 , p + £ n (p, 9, yi ) > -B 2 }, 


k -1 

H L(f? fc>n )" 1 L(f? fc>n + a(R k , n ){p + Up, 0, Vi))) < C 2 (p V l)^” 1 ). 

2= 1 

Since a is eventually non-increasing, 


fc-i 

expj U [R ktn + a(R k>n ) (p + (p, 0, tp)) ) + <tj)(Rk , n )} 

/e-1 

< exp{-(p + £ n (p,0,?p))} . 

2=1 


Now, we have 

n fc P{/i n (44) = 1, ll-Xjll - R n > —a(Rk, n )B 2 , i = l,...,k} 

k _1 

< 4CC 2 f [ [ (pV JJexp{-(p + £ n (p,0,y;))) 

x 1 (p + £ n (p, 9, yi) > -B 2 ) h( 0, y)dpJ ( 6)d6dy . 
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If we can show that 

(6.15) expj-(/9 + £ n (/9,0,yi)) j l(p + £ n (p, 9, y t ) > -B 2 ) — > 0 , n —y oo 

for every p > —B 2 , 9 G S d 1 , and |y, || < M, i = 1 , d, then the dominated 
convergence theorem finishes the proof. First, in the case of (6,yi) < 0, 

P + £n(p, 9, Vi) -oo, rt-> oo 


and hence, as n —>• oo, 

exp{-(p + &i(p,0,yi))} l(p + Zn(p,6,yi) > ~B 2 ) 

< e B2 1 {p + €n(p,0,yi) > -B 2 ) 0 . 

Second, if (9,yi) > 0, then exp{ — [p + £ n (p, 9, j/j))} —)• 0. So in either case, 
(6.15) is established. □ 
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